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Abstract 

We review recent results for the properties of doped antiferromagnets, obtained by the nu- 
merical analysis of the planar t-J model using the novel finite-temperature Lanczos method for 
small correlated systems. First we shortly summarize our present understanding of anomalous 
normal-state properties of cuprates, and present the electronic phase diagram, phenomenologi- 
cal scenarios and models proposed in this connection. The numerical method is then described 
in more detail. Following sections are devoted to various static and dynamical properties of the 
t-J model. Among thermodynamic properties the chemical potential, entropy and the specific 
heat are evaluated. Discussing electrical properties the emphasis is on the optical conductivity 
and the d.c. resistivity. Magnetic properties involve the static and dynamical spin structure 
factor, as measured via the susceptibility measurements, the NMR relaxation and the neutron 
scattering, as well as the orbital current contribution. Follows the analysis of electron spectral 
functions, being studied in photoemission experiments. Finally we discuss density fluctuations, 
the electronic Raman scattering and the thermoelectric power. Whenever feasible a comparison 
with experimental data is performed. A general conclusion is that the t-J model captures well 
the main features of anomalous normal-state properties of cuprates, for a number of quantities 
the agreement is even a quantitative one. It is shown that several dynamical quantities exhibit 
at intermediate doping a novel universal behaviour consistent with a marginal Fermi-liquid con- 
cept, which seems to emerge from a large degeneracy of states and a frustration induced by 
doping the antiferromagnet. 
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1 Introduction 



The discovery of the high-temperature superconductivity (SC) in copper-oxide based compounds - 
cuprates (Bednorz and Miiller 1986) revived the interest in materials containing transition elements. 
The main feature of these materials is the crucial role of electron-electron interactions. This can 
result in electronic properties, very unusual when compared to the behaviour of conventional metals. 
Although the unconventional SC is clearly the most puzzling phenomenon in cuprates, we focus in 
this review predominantly on the analysis of more or less anomalous properties of the normal state 
which deviate essentially from the standard understanding of electrons in metals and still present 
one of the major theoretical challenges in the solid state physics. 

The cuprate superconductors have a very anisotropic structure, where the common building 
blocks are layers, in cuprates formed by combining one of the three possible structural elements 
containing Cu and O, as shown in Fig. |l.l[b. The CuO layered structures are stacked in the crystal, 
separated by various intercalant layers in different cuprates. In spite of vast differences in the 
structure of unit cells, electronic properties of the whole family of cuprates are quite universal. 
This can be explained by the predominant role of generic CuC-2 planes, Fig. |l.l| a, where conducting 
electrons reside. The electronic coupling between Cu02 planes is very weak, resulting in a huge ratio 
of the in-plane resistivity to the perpendicular resistivity, in the anisotropy of the SC coherence 
length etc. Intermediate layers serve mainly as a charge reservoir for the planes. Consequently 
properties of cuprates are quite well classified according to the doping level of the reference Cu02 
electronic structure. 




Figure 1.1: (a) Schematic structure of copper-oxide planes and (b) three possible building blocks 
of the planes, after Fulde (1991). 

As now well established, the reference cuprate compounds, as La2Cu04 and YBa2Cu306, are 
Mott (charge-transfer) insulators due to strong correlations which induce in a half-filled band 
a charge gap ~ 2 eV. The spin degrees of planar Cu02 electrons can be well mapped on the 
properties of a planar antiferromagnetic (AFM) 5 = 1/2 Heisenberg model. The AFM ground 
state of undoped materials, emerging from strong correlations, has been quite early recognized as 
a crucial starting point for theoretical considerations (Anderson 1987) of doped materials, being 
strange metals in the normal state and exhibiting high transition temperature T c to the SC state. 

There are by now numerous indications that essential features of electronic properties of the 
doped AFM, as realized in cuprates, are well represented by prototype single-band models of 
correlated electrons, as the Hubbard model and the t-J model (Rice 1995). In spite of their 
apparent simplicity both models are notoriously difficult to treat analytically, in particular in 
the most interesting regime of strong correlations. The lack of analytical tools for correlated 



4 



electrons (for a general introduction see Fulde 1991) has increased the efforts towards numerical 
approaches (Dagotto 1994), which predominantly can be divided into two categories: the quantum 
Monte Carlo (QMC) methods and the exact diagonalization (ED) methods. The t-J model, which 
incorporates the strong correlation requirement explicitly, is more adapted to the ED approach. 
So far most calculations were performed for the ground state (g.s.) at T = 0, where the standard 
Lanczos technique (Lanczos 1950) offers an efficient ED analysis of small systems of reasonable 
sizes. Recently, present authors (Jaklic and Prelovsek 1994a) introduced a novel numerical method, 
combining the Lanczos method with a random sampling, which allows for an analogous treatment of 
statics and dynamics of many-body quantum models at T > 0. The latter method, further referred 
to as the finite-temperature Lanczos method (FTLM), and results for the t-J model obtained using 
this method, are the main subject of this review. 

The final goal is to understand properties of doped AFM in general, and of high-T c cuprates 
in particular. In the absence of reliable analytical methods and results, numerical calculations 
can help to answer several crucial questions. Are strong correlations, as incorporated in prototype 
models, enough to account for anomalous normal-state properties of the strange metal ? Which 
are the relevant energy and temperature scales in doped AFM, as represented by the t-J model, 
and in which properties do they show up ? Which is the unifying phenomenological description of 
the normal state ? Is the t-J model sufficient, or which ingredients should be added to account 
qualitatively and quantitatively for observed properties ? Can we learn something macroscopically 
meaningful from the studies of small systems, and why ? 

We note that at present the SC seems to be beyond the reach of numerical approaches including 
the FTLM, hence we do not investigate here in more detail the possible existence of SC and its 
origin in model systems. 



2 Cuprates as doped antiferromagnets 

2.1 Electronic phase diagram and properties of the normal state 

Reference cuprate compounds are AFM insulators and are so far best understood. Properties of 
various other layered cuprates can be interpreted in terms of doping the reference material, where 
(mainly) holes are introduced into Cu02 planes. One of the major conceptual achievements, which 
emerged from careful experimental investigations of high-quality materials in the last decade, has 
been the realization of quite universal electronic phase diagram (Hwang et al. 1994, Batlogg et 
al. 1994, Batlogg 1997), revealing characteristic temperature scales as they develop as the function 
of hole doping. It is at present quite common to classify materials with respect to doping as 
underdoped, optimally doped and overdoped, and the corresponding phase diagram is shown in 



Fig. 2.1 



2.1.1 Optimum doping regime 

Experimentally optimum doping is chosen to correspond to materials with the highest T c within 
the given class of chemically and structurely related compounds. It has been realized soon after 
the discovery of high-T c SC that also normal state properties at T > T c of the optimally doped 
materials are very anomalous, but at the same time also most universal. The prominent feature 
is the resistivity law p oc T (Takagi et al. 1992), valid essentially in the whole measurable range 
T > T c and clearly contradicting the normal Landau- Fermi- liquid (LFL) behaviour p oc T 2 . Related 
is the observation that the dynamical conductivity ct(lu) does not fall off for larger ui according to 
the Drude form a oc 1 /uj 2 , but rather as a oc 1/u (Tanner and Timusk 1992). 
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Hole Doping 



Figure 2.1: Schematic electronic phase diagram of cuprates, after Batlogg (1997). 



The clearest evidence for an anomalous spin dynamics comes from the NMR and NQR relaxation 
(Slichter 1994), where the relaxation rate 1/Ti on planar 63 Cu in optimally doped La2_ x Sr x Cu04 
(LSCO) with x opt ~ 0.15 is nearly T-independent (as well as nearly doping independent) in contrast 
to the usual Korringa law for metals \ jT\ oc T. The qualitative difference, i.e. a large enhancement 
of low-frequency spin fluctuations at low T, can be related to the persistence of short range AFM 
fluctuations, even at the optimum doping. The support for this comes also from neutron scattering 
experiments, where a substantial AFM correlation length £ has been measured in the same class of 
materials, with £ ~ 3.8A/^/x (Birgeneau et al. 1988). 

On the other hand, several properties at the optimum doping seem to be close to the normal LFL 
picture. The angle-resolved photoemission spectroscopy (ARPES) on cuprates (Shen and Dessau 
1995), most reliable for BiSrCaCuO (BISCCO) compounds, shows electronic excitations consistent 
with a large Fermi surface (FS) and with the conserved FS volume (Luttinger theorem), although 
spectral shapes are quite distinct from a simple LFL picture. The specific heat (Loram et al. 1993) 
follows in the normal state roughly the LFL behaviour C v = 7T with 7 not very far from the 
free-fermion value, consistent with a nearly T-independent uniform susceptibility xo> whereby the 
Wilson ratio is close to the free-fermion one (Loram et al. 1996). It is also the unifying characteristic 
of the optimum doping that properties do not reveal above T c any additional characteristic T scale. 



2.1.2 Overdoped and underdoped regime 

In overdoped materials T c is decreasing and finally vanishing with increased doping. At the same 
time the electronic properties are getting closer to the usual metallic behaviour consistent with the 
LFL scenario. E.g., the resistivity behaviour moves towards the normal FLF form p oc T 2 (Takagi 
et al. 1992), spectral shapes of electronic excitations, as revealed by ARPES, become sharper near 
the FS (Marshall et al. 1996) etc. These facts can be put together with a decreasing intensity 
of AFM fluctuations. It is thus plausible that we are dealing in the overdoped regime with the 
crossover to the normal LFL, however this crossover is not a trivial one and so far also not well 
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understood either. 

The most evident progress in the investigations of normal-state properties has been made in 
last few years for the underdoped cuprates. In contrast to the optimum doping, experiments reveal 
in this regime additional characteristic temperatures T > T c (Batlogg et al 1994), which show 



up as the crossovers where particular properties qualitatively change. As summarized in Fig. |2J. 
there seems to be an indication for two distinct crossovers. The existence of both as well as their 
distinction is still widely debated, nevertheless we will refer to them as the AFM crossover scale 
T* and the pseudogap scale T sg for the lower one (Batlogg 1997). 

The T* scale (Batlogg et al. 1994) shows up most clearly as the maximum of the spin suscepti- 
bility xo(T = T*) (Torrance et al. 1989). The in-plane resistivity p{T) is linear p oc T for T > T* 
and decreases more steeply for T < T*. Characteristic is also the anomalous T-dependence of the 
Hall constant R H {T) for T < T* (Ong 1990, Hwang et al. 1994). The latter is evidently hole-like 
in the classical sense Rh(T > T c ) oc 1/x, which is also not properly understood theoretically. It 
seems plausible that the T* crossover is related to the onset of short-range AFM correlations for 
T < T* , since in the undoped AFM T* corresponds just to a well understood maximum due to a 
gradual transition from a disordered paramagnet to the one with short-range AFM correlations. 

The crossover T S g has been first identified in connection with the decrease of the NMR relaxation 
1/T\ for T < T sg (Takigawa et al. 1991, Slichter 1994), indicating the reduction of low-energy spin 
excitations interpreted as the opening of the spin pseudogap in underdoped materials. Most striking 
evidence for an additional energy scale in underdoped cuprates is the observation of the leading- 
edge shift in ARPES measurements at T > T c (Marshall et al 1996), indicating features of the 
d-w&ve SC gap persisting within the normal phase. It should be pointed out, however, that the 
designation of crossover features is at present still controversial. In particular it is not evident 
whether we are dealing with two or more essentially different energy scales. 

2.2 Phenomenology of the normal state 

Properties of the normal LFL follow from the one-to-one correspondence of low-energy excitations 
of the interacting fermion system to that of a free-fermion gas. The prerequisite is that the volume 
of the FS is conserved (Luttinger 1960). Essential are the well defined quasiparticles (QP) with 
the vanishing damping at the FS, with r oc (E — Ep) 2 . Consequences are the linear specific heat 
Cy = 7T, nearly T-independent static and dynamical spin susceptibilities, the Korringa law for 
NMR relaxation 1 jT\ oc T, the resistivity p(T) oc T 2 etc. Experimental facts on cuprates contradict 
the usual LFL picture. Several more or less elaborated scenarios have been proposed to capture 
main anomalous features. 

Focusing on the importance of AFM spin fluctuations, the concept of a nearly AFM Fermi 
liquid (NAFL) has been elaborated (see e.g. Monthoux and Pines 1994). Here one assumes that at 
low-frequencies x"(<f, w) oc u at all q, as expected in a LFL, but with strongly enhanced fluctuations 
near the AFM wavevector Q = (w, 7r), corresponding to the critical slowing down in the proximity 
of a phase transition to the AFM-ordered state. The following form has been proposed, which can 
be derived also via the self-consistent paramagnon theory (Moriya et al. 1990), 

x(M = ^ • (2-1) 

l + ?\q-Q\*-iu/u SF 

In the proposal by Millis et al. (1990), originally devoted to the interpretation of the NMR and 
NQR relaxation, the main T-dependence is expected to arise from the AFM correlation length £, 
which is assumed to show a critical behaviour as T is decreased, i.e. £ 2 oc T X /(T + T x ). In order 
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to explain the anomalous NQR relaxation \/T\(T) ~ const (Imai et al. 1993), contradicting the 
Korringa law, strong T-dependence of £ is essential with T x ~ 100 K, as well as low oj$f and large 
Xq for T ~ T c (Monthoux and Pines 1994). The form of spin fluctuations is the basis for further 
investigations of the charge dynamics, strongly coupled to spin degrees of freedom. The calculated 
response functions such as the conductivity appear also anomalous, e.g. the resistivity is close to a 
linear law p oc T. 

There are several proposals, analogous to the NAFL in the basic idea of the proximity to a critical 
point, enhancing fluctuations in the optimum doping regime. One proposal invokes the quantum 
critical scaling of the spin dynamics, established in the nonlinear sigma model (Chakravarty et al. 
1989), induced by doping the AFM (Sokol and Pines 1993). Another scenario relates the quantum 
critical point to a charge-density- wave instability (Castellani et al. 1995). 

An alternative interpretation of experimental facts has been provided by the concept of a 
marginal Fermi liquid (MFL) (Varma et al. 1989). The hypothesis is that there exist excitations, 
contributing both to the charge and spin response, which show in a broad range of wavevectors q 
anomalous susceptibilities of the form 

///—* \ f C uj/T for \u\ < T, 
X ~ < . . (2.2) 

[ G sgn lo for > T 



Due to the scattering on bosonic excitations with the spectrum ( |2.2D the single-particle self-energy 
Yi(k,uj) is also anomalous. Assuming its ^-independence in a broad range, it has been postulated 
using phenomenological arguments (Littlewood and Varma 1991), 



S(w) ~ 7rA 



X 

2lu In ix 



x = max(\u\,nT), (2.3) 



where lo c is a high-frequency cutoff. Hence the QP lifetime 1/r oc —T,"(uj) is anomalous, i.e. 
1/t(cj) oc x. It should be however mentioned that the Ansatz ( |2.3| ) is not unique and also modified 
forms have been invoked, e.g. in the analysis of the optical conductivity (El Azrak et al. 1994, 
Baraduc et al. 1995) better fit has been obtained with 

-E"(w) ~vrA(M +ttT). (2.4) 

While the FS should remain well defined with the volume equal to that of free fermions, the 
corresponding QP weight Z at the FS, given by 

Z- 1 = l-dE'(u)/du\ w=0 , (2.5) 

vanishes on the FS in analogy to the case of a one-dimensional Luttinger liquid (Haldane 1981). The 
MFL concept accounts for several remarkable properties at the optimum doping, as the anomalous 
resistivity p(T) oc T, the optical conductivity o~(uj) oc the NMR and the NQR relaxation rate 
\/T\(T) ~ const. Note that the only low-energy scale within the MFL scenario, equations ( |2.2| ) - 
(|2.4| ), is given by T. Although there are certain similarities to the NAFL and other critical-point 
scenarios, the essential difference of the MFL concept is a non-critical q dependence. Hence the 
critical behaviour within the MFL is rather a local one. The attempts to derive MFL behaviour 
from a microscopic model have not been successful so far. 

2.3 Models of correlated electrons in cuprates 

A similarity of electronic properties in a wide class of different cuprates serves as a strong indication 
that the appropriate microscopic model should be quite universal and must in first place describe 
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the electrons restricted to CuC>2 orbitals within a single plane. There has been quite an extensive 
effort put into finding a proper model, and at present there seems to be a wide consensus on its 
main features. This should be contrasted to various other materials with interacting electrons, e.g. 
heavy fermions and ID conductors, where microscopic models are much less known. 

Since the physics of electrons in Cu02 planes is governed by Cu 2>d x 2_ y 2 orbitals and O 2p X;V 
orbitals (see the structure in Fig. |1.1| ), quite a complete model seems to be the three-band Hubbard 
model (Emery 1987), describing fermion carriers (holes) added to closed 3d 10 and 2p 6 shells. Pa- 
rameters are the Cu-0 hopping t pc i, the direct 0-0 hopping t pp , the on-site energies Ed, £ P , and the 
corresponding Coulomb repulsions Ud, U p on Cu and O sites, respectively. Parameters correspond 
to the charge-transfer regime with A = e p — Ed > t p d and Ud > A. The reference (insulator) 
material contains one fermion/cell, entering predominantly the d orbitals. Due to Ud S> t p d, a large 
charge gap ~ 2eV opens at half filling, while spin degrees can be mapped on the isotropic S = 1/2 
Heisenberg model, first proposed in connection with cuprates by Anderson (1987). Holes added by 
doping enter the singlets, introduced by Zhang and Rice (1988), which can be in a fermion model 
treated as empty sites (holes). Such a reduction, confirmed also with other analytical approaches 
(Zaanen and Oles 1988, Ramsak and Prelovsek 1989) and cluster methods (Hybertsen et al. 1990), 
leads to a single-band t-J model (Rice 1995), 

H = - t J2(c} s Cis + H.c.) + J ^(5, • Sj - \n % n 3 ), (2.6) 

{ij)s (ij) 

describing fermions in a tight-binding band with the hopping parameter t oc t^ d j 'A. Here Si = 

(V^)X)ss' Hs^ss'Cis 1 are the local spin operators interacting with the exchange parameter J oc 
tp d /UdA < t. Due to the strong on-site repulsion states with doubly occupied sites are explicitly 
forbidden and we are dealing with projected fermion operators C{ s = Ci s {\ — nj._ s ). 

By explicitly projecting out states with doubly occupied sites, the t-J model only allows for 
charge fluctuations in terms of a hole motion, while at half-filling it is equivalent to the S = 1/2 
Heisenberg model. The t-J model is expected to capture the essential low-energy physics of doped 
AFM as well of cuprates in the whole regime of dopings. The challenging regime of the model is 
the one of strong correlations with J < t. 

The t-J model is the simplest model which describes the interplay of magnetism and itinerant 
metallic properties of cuprates. A more rigorous reduction of the three-band model (Zaanen and 
Oles 1988, Ramsak and Prelovsek 1989) leads to additional terms, which could be as well repre- 
sented within a reduced space without doubly occupied sites. Among possible generalizations most 
attention has been recently devoted to the addition of the next-nearest-neighbour (n.n.n.) hopping 
{t') terms, emerging from the t pp hopping in the three-band model, 

H v = -t' (c] s c is + R.c), (2.7) 

««»« 

representing the hopping along the diagonal in Fig. |l.l| a. Analogous is the t" term for the n.n.n. 
hopping along each axis. It seems necessary to include such a term to account for the QP dispersion 
found in ARPES experiments on undoped material, such as Sr2Cu02Cl2 (Wells et al. 1995). In 
spite of their apparent smallness, t' and t" terms could lead to relevant corrections since they allow 
a free propagation of fermions even in an AFM (Neel) spin background. 

The t-J model, as relevant for cuprates, should be considered on a planar square lattice. Pa- 
rameters are rather well known (Rice 1995). J is measured in the undoped AFM via the inelastic 
neutron scattering and the magnon dispersion, leading to J ~ 0.13 eV. The hopping parameter t 
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is not accessible directly, but cluster calculations (Hybertsen el al. 1990) and other considerations 
(Rice 1995) allow only for a narrow range of t values. In our calculations we shall furtheron use (if 
not declared differently) t = 0.4 eV and J/t = 0.3. A possible range for t' (Hybertsen et al. 1990) 
is more controversial, while in numerical studies (Tohyama and Maekawa 1994, Nazarenko et al. 
1995) values -0.35 < t'/t < -0.2 are used. 



Another prototype model for strongly correlated electrons is the traditional Hubbard model 
(Hubbard 1963), 



High-energy excitations of the Hubbard model could be different from those of the charge-transfer 
regime in the three-band model, still it is expected that the low-energy properties map well on those 
of the t-J model provided that U S> t. In the following we shall mainly consider T > properties 
of the t-J model. It should be however noted that prior to the introduction of the FTLM method 
most calculations of T > properties have been performed for the Hubbard model by applying 
QMC methods (Dagotto 1994). 

3 Finite- Temperature Lanczos Method 

This chapter is devoted to the description of the FTLM which we developed (Jaklic and Prelovsek 
1994a) for studying correlated systems at T > and which is used to obtain results described 
furtheron. The goal was to calculate T > properties in small model systems and to find a 
method, comparable in efficiency to g.s. calculations employing ED methods, used in the past 
decade extensively in the study of correlated systems (Dagotto 1994). 

Here we should stress that the advantage of T > calculations is twofold. It is evident that we 
are interested in static and dynamical properties at nonzero T, in particular in their T-variation. 
On the other hand, the use of finite but small T > represents the proper approach to more reliable 
g.s. calculations in small systems. Namely, it is well known that g.s. ED results, in particular for 
dynamical quantities, are strongly influenced by finite size artifacts. At T > the latter effects can 
to large extent average out, leading to more macroscopic-like results. Still the understanding of 
remaining finite-size restrictions is important for the proper application of the method, as will be 
described in Sec. 3.7. In Sec. 3.8. we put our approach in perspective with other methods yielding 
T > results for models of correlated electrons. These includes mainly various QMC methods and 
the high-T expansion (HTE) technique. 

3.1 Lanczos algorithm and matrix elements 

The scarcity of well-controlled analytical approaches to models of strongly correlated electrons has 
stimulated the development of computational methods. Conceptually the simplest is the ED method 
of small systems. In models of correlated electrons, however, one is dealing with the dimension of 
the basis (number of basis states) which grows exponentially with the size of the system. In the 
Hubbard model there are 4 basis states for each lattice site, therefore the number of basis states 
in the iV-site system is N st oc 4 N . In the t — J model N st still grows as oc 3^. In the ED of such 
systems one is therefore representing operators with matrices N s t x N s t, which become large already 
for very modest values of N. The helpful circumstance is that for most interesting operators and 
lattice models only a small proportion of matrix elements is nonzero within the local basis. Then, 
the operators can be represented by sparse matrices with N st rows and at most f(N) nonzero 
elements in each row. In this way memory requirements are relaxed and matrices up to N st ~ 10 7 




(2.8) 
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are considered in recent applications. Finding eigenvalues and eigenvectors of such large matrices is 
not possible with standard algorithms performing the full diagonalization. One must instead resort 
to the power algorithms (see Parlett 1980), among which the Lanczos algorithm (Lanczos 1950) is 
one of the most widely known. 

The Lanczos algorithm starts with a normalized random vector |0q) in the vector space in which 
the Hamiltonian operator H is defined. H is applied to |<fo) an d the resulting vector is split in 
components parallel to |(/>o), and orthogonal to it, respectively, 

#|0 o )=ao|0 Q ) + &i|0i). (3.1) 

Since H is Hermitian, ao = (^ol-^l^o) is real, while the phase of can be chosen so that b\ is 
also real. In the next step H is applied to \4>\), 

i?|0i> = &i]0 o ) + + 6 2 |02>, (3.2) 

where |</>2) is orthogonal to |0o) and \4>i). It follows also b[ = {cpo\H\4>i) = b\. Proceeding with the 
iteration one gets in % steps 

E\<\>i) = 6 l -|0i_i) + a i |0i)+6 m |0 i+ i), l<i<M. (3.3) 

By stopping the iteration at i = M and putting the last coefficient &Af+i = 0, the Hamiltonian can 
be represented in the basis of orthogonal Lanczos functions \4>i) as the tridiagonal matrix Hm with 
diagonal elements ai with i = . . . M, and offdiagonal ones bi with i = 1 . . . M. Such a matrix is 
easily diagonalized using standard numerical routines to obtain approximate eigenvalues €j and the 
corresponding orthonormal eigenvectors 

M 

hk-> = l>ji|&>, j = 0...M. (3.4) 

It is important to realize that are (in general) not exact eigenfunctions of H, but show a 
remainder 

H\ipj) - ej\if)j) = b M+1 v jM \(t>M+i}- (3.5) 
On the other hand it is evident from the diagonalization of Hm that matrix elements 

(ip i \H\?p j ) = e j S ij , i,j = 0...M (3.6) 

are exact, without restriction to the subspace Lm- 

If in the equation (3.3) &a/+i = 0, we have found a (M + l)-dimensional eigenspace where Hm 



is already an exact representation of H. This inevitably happens when M = N s t — 1, but for 
M < — 1 it can only occur if the starting vector is orthogonal to some invariant subspace of H. 
This should not be the case if the input vector |^>q) is random, without any hidden symmetries. 

The number of operations needed to perform M Lanczos iterations scales as MN s t- Numerically 
the Lanczos procedure is subject to roundoff errors, introduced by the finite-precision arithmetics. 
This problem usually only becomes severe at larger M > 100 (more than needed to get accurate 
g.s. IV'o)) an d is seen in the loss of the orthogonality of vectors \<pi)- It can be remedied by 
successive reorthogonalization (and normalization) of new states \4>i), plagued with errors, with 
respect to previous ones. However this procedure requires ~ M 2 N s t operations, and can become 
computationally more demanding than Lanczos iterations alone. This effect prevents one to use 
the Lanczos method e.g. to tridiagonalize large matrices. 
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The identity fl3.6| ) already shows the usefulness of the Lanczos method for the calculation of 
particular matrix elements. As an aid in a further discussion of the Lanczos method we consider 
the calculation of a matrix element 



W, 



ki 



(n\H k BH l A\n) 



(3.7) 



where \n) is an arbitrary normalized vector, and A, B are general operators. One can calculate this 
expression exactly by performing two Lanczos procedures with M = max(fc, I) steps. The first one, 
starting with the vector \(J)q) = \n), produces the subspace Lm = {\<l>j), 3 = . . . M} along with 
approximate eigenvectors \ipj) and eigenvalues €j. The second Lanczos procedure is started with 
the normalized vector 

\4> ) = A\<j) Q )/^MMA^), (3.8) 

and results in the subspace Lm = {|</>j), j = . . . M} with approximate \tpj) and ij. We can now 
define projectors 



8=0 

which for m = M can also be expressed as 

M 



i=0 



M 



i=0 



i=0 



By taking into account definitions ( |3.9| ), ( |3.10 ) we show that 



HP„ 



Pm+ 1 PtPm 



PAlHP m , 



m < M. 



(3.9) 



(3.10) 



(3.H) 



6q) = Pol^o) and A\n) oc |0o) = Po\(^o), by successive use of the first 



Since in addition |n) : 
equality in ( 3.11 ) we arrive at 

W kl = {MPoHPtH . . . HP k BPiH . . . P x HhA\<fo). 



(3.12) 



Using the second equality in the equation ( 3.11 ) and identities -Pol^o) = Pm\<Po), PoA\4>q) 
PmA\(J)q) we can rewrite Wja as 



Wi 



ki 



b \P M HP M H . . . HP M BP M H . . . P M HP M A\<j) Q 



(3.13) 



We note that the necessary condition for the equation ( 3.13 ) is M > k,l. We finally expand the 
projectors according to expressions ( |3.10D and take into account the orthonormality relation ( p.6[ ) 
for matrix elements, and get 

M M M M 

W kl = ^...^^...^(0o|^ o )(V'ioWi 1 )---(^-J^iJ 

«o=0 ifc=Ojo=0 ji=0 

xfak\ B $h)$ji\ H \'>i>3t-i) ■ ■ ■ $h\H$j )$j \A\<h) = 

M M 

We have thus expressed the desired quantity in terms of the Lanczos (approximate) eigenvectors, 
and eigenvalues alone. 
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3.2 Dynamical response in the ground state 

Within the Lanczos algorithm the extreme (smallest and largest) eigenvalues ej, along with their 
corresponding are rapidly converging to exact eigenvalues Ei and eigenvectors It is quite 
characteristic that usually (for nondegenerate states) M = 30 — 60 <C N s t is sufficient to achieve 
the convergence to the machine precision of the g.s. energy E$ and the wavefunction l^o)! horn 
which various static and dynamical correlation functions at T = can be evaluated. 

After l^o) is obtained, the g.s. dynamic correlation functions can be calculated within the same 
framework (Haydock et al. 1972). Let us consider the autocorrelation function 

C(t) = -i{V Q \A\t)A\$ Q ) = -i(^ \A^e^ E °- H ^A\^ ) (3.15) 

with the transform, 

/•OO 1 

C(u) = / dte^ +t C(t) = <^o|At_ r _— -A\V ). (3.16) 

where uj + = uj + ie, e > 0. To calculate C(u), one has to run the second Lanczos procedure 
starting with the normalized function \4>o), equation (3.8). The matrix for H in the new basis Lm, 
with elements {(j>i\H\(j>j) = [Hj\f]ij, is again a tridiagonal one with Oj and hi elements, respectively. 
Terminating the Lanczos procedure at given M, one can evaluate the C{u) as a resolvent of the 
Hm matrix which can be expressed in the continued- fraction form (Haydock et al. 1972), 

CH = gg^W p , (3.17) 

lo+ + E - a T 

^ H 

uj+ +E - ai 



lo+ + E - a 2 



terminating with &m+i = 0, although other termination functions have also been employed. 
The spectral function C(uS) = — (l/7r)ImC(w) is characterized by frequency moments, 



/oo 
Jc(u)du = (* \Ai(H - ^0)^1^0), (3-18) 
-00 

which are particular cases of the expression Q3.7| ) for B = A', k = 0, and \n) = Using the 

equation ( |3.14| ) we can express \i\ for I < M in terms of Lanczos quantities alone 

M 

m = ^(* |At|^)(^|A|*o>(€i - E ) 1 . (3.19) 

j=0 

Hence moments are exact for given |^o) ( as would be also for any other starting \n)) provided 
I < M. The corresponding approximation for C(t), equation ( 3. 15| ) , within the restricted set of 
eigenfunctions \ipj),j = . . . M, can be written at given M (assuming &a/+i = 0) as 

M 

C(t) = -^(^ol^l^-X^-IA^o)^-^)*. (3.20) 
3=0 

Note that such C{t) expanded as a series in t (short-time expansion) has exact M terms, since the 
coefficients are just moments fii, equation ( p.!9[ ). 
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As a practical matter we note that |(V5j|yl|^ / o)| 2 = ^Jo> nence no matrix elements need to be 
evaluated within this approach. In contrast to the continued fraction fl3.17| ), the expression ( 3.20 ) 



allows also the treatment of more general correlation functions (B(t)A), with B ^ A' . In this case 
the matrix elements {*$>o\B\vpj) have to be evaluated explicitly. From the above one sees that the 
Lanczos method is very convenient to calculate the frequency moments as well as the dynamical 
C(lo). Certain ideas presented above will be used to construct the algorithm for T > 0, discussed 
in next subsections. 

3.3 High-temperature expansion 

The novel method for T > is based on the application of the Lanczos iteration, reproducing 
correctly high-T and large-u; series. The method is then combined with the reduction of the full 
thermodynamic trace to the random sampling. We present these ingredients in the following. 
We first consider the expectation value of the operator A in the canonical ensemble 

N st , N st 



(A) = Y,(n\e' /3H A\n) J ^(n|e^ H |n), (3.21) 



n=l ' n=l 



where (3 = 1/ksT. A straightforward calculation of (A) requires the knowledge of all eigenstates 
|\I/ n } and corresponding energies E n , obtained by the full diagonalization of H, 



N at , N s - 



(A) = Y j e- f3En (*n\A\-Z n ) / J>-^, (3.22) 



n=l ' n=l 



computationally accessible only for N s t ~ 5000. Instead let us perform the HTE of the exponential 



e-P H , 

Nst oo , R -. k 
n=l k=0 

N st OO , R , k 

n=l fc=0 

Terms in the expansion (n\H k A\ n) can be calculated exactly using the Lanczos procedure with 
M > k steps and with \4>q) = \n) as a starting function, since this is a special case of the expression 
Using the relation ( |3.14j ) with I = and B = 1, we get 

M 

(n\H k A\n) = ^(n|^)(Vr|A|n)(e?) fe . (3.24) 
i=0 



Working in a restricted basis k < M, we can insert the expression (3.24) into su ms (|3.23| ), extending 



them to k > M. The final result can be expressed in analogy to the equation (|3.2C| ) as 

N at M 

n=l i=0 

Nst M 

z w EE e_w ( n W)wi ft )- ( 3 - 25 ) 

n=l i=0 
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and the error of the approximation is of the order of f5 M+1 . 



Evidently, within a finite system the expression (|3,25|) , expanded as a series in (3, reproduces 
exactly the HTE series to the order M. In addition, in contrast to the usual HTE, it becomes 
(remains) exact also for T — > 0. Let us assume for simplicity that the g.s. |^o) is nondegenerate. 
For initial states |n) not orthogonal to l^o)) already at modest M ~ 50 the lowest function \ipo) 
converges to l^o)- We thus have for j3 — > do, 



N si , N ax 

n=l 1 n=l 

<*o|A|* )/<*o|*o), (3-26) 



(A) = £(n|tf ><¥o|4n>/£(n|*o>(*o| 



where we have taken into account the completeness of the set \n). Obtained result is just the usual 
g.s. expectation value of an operator. 

3.4 Large-frequency expansion at T > 

In order to calculate dynamical quantities, the HTE must be supplemented by the high-frequency 
(short-time) expansion analogous to the one used at T = in deriving the equation ( |3.20| ) from 
( |3.19| ). The goal is to calculate the dynamical correlation function at T > 0, 



(B(t)A) = Tr e- pil e int Be- int A /Tire - *". (3.27) 
Expressing the trace explicitly and expanding the exponentials, we get 

(B(t)A) =Z-^f:f: ti + tpl (n\H k BH l A\n). (3.28) 

n=l fe=0 1=0 

Expansion coefficients in equation (|3.28j ) can be again obtained via the Lanczos method, as discussed 
in Sec. 3.1. Performing two Lanczos iterations with M steps, started with normalized \4>q) = \n) 



and \<Pq) oc A\n), respectively, we calculate coefficients Wm following the equation (|3.14[) , while Z 
is approximated by the static expression ( 3.25| ). Extending and resumming series in k and I into 
exponentials, we get 

N at M M 

(B(t)A) « r^^^e^e^-^^ICWlBI^Xtl^k). (3-29) 

71=1 4 = j=o 

We check again a nontrivial T = limit of the above expression. If |n) are not orthogonal to 
the g.s. |^o)> then for large enough M the lowest-lying state converges to ~ Eq and \^)q) ~ l^o), 
respectively. In this case we have in analogy to the equation ( |3.26| ) 

Nst M I 

(B(t)A) « ^E^^^l^l^) A' l^l^o>(<^ol^o> / <*o|*o>. (3-30) 

71=1 j = o ' 

Generally larger M 3> 100 are needed in order that relevant higher-lying states |^") and become 
independent of |n). Only in such a limit we recover strictly the g.s. result, corresponding for 
A^ — > B to the equation ( 3. 16| ) ■ Note however that similar restrictions apply to the continued 



fraction expansion (|3.17|) which reproduces correctly moments {i\ (|3.18| ) up to I < M, but not 
necessarily the details (e.g. positions and weights of peaks) of the C{oS) spectrum. 
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3.5 Random sampling 



The computation of static quantities ( |3.25| ) and dynamical ones ( p. 29 ) still involves the summation 
over the complete set of N st states |n), which is not feasible in practice. To obtain a useful method, 
one further approximation must be made which replaces the full summation by a partial one over a 
much smaller set of random states (Imada and Takahashi 1986). Such an approximation analogous 
to Monte Carlo methods is of course hard to justify rigorously, nevertheless we can estimate the 
errors involved. 



We consider the expectation value {A) at T > 0, as defined by the expression ( 3.21 ). Instead 
of the whole sum in equation ( 3.21 ) we first evaluate only one element with respect to a random 
state |r), which is a linear combination of basis states 



N s 



n=l 



n) 



(3.31) 



i.e. P rn are assumed to be distributed randomly. Let us discuss then the random quantity 

A r = {r\e-P H A\r)J{r\e-P H \r) = 

N st , N st 

= P*mPrm(n\e- pH A\m) / ff n p rm (n\e-P H \m). (3.32) 



n,m=l 



n,m=l 



We choose for convenience that here basis states |n) correspond to the eigenstates of H. We first 
assume in addition that [H, A] = 0, to diagonalize simultaneously both A and H. Then we have 



A r = Y\Pm\ 2 (n\e-f SH A\n) / £ \^ n \ 2 {n\e^ H 



n) 



(3.33) 



n=l 



n=l 



We can express |/? r n| 2 = 1/^Vst + <5 rn , where random deviations 5 Tn are not correlated with matrix 
elements (n\e~^ H \n) = Z n and {n\e~^ H A\n) = Z n A n . It is then easy to see that A r is close to (^4), 
and the statistical deviation is related to the effective number of terms Z in the thermodynamic 
sum, i.e. 

Ar = {A) + 0(1/ V^), Z = e^ Eo Y Z n . (3.34) 

n 

Note that for T — > oo we have Z — * N s t and therefore at large N s t a close estimate of the average 
( |3.34| ) can be obtained from a single random state (Imada and Takahashi 1988, Silver and Roder 
1994). On the other hand, at finite T < oo the statistical error of A r increases with decreasing 
Z. Still, strictly at T = and for a nondegenerate g.s. we obtain from equation ( |3.33| ) again the 
correct result. 

In the FTLM we replace the full summation in the expression ( 3.2 1] ) with a restricted one over 
several random vectors |r), r = 1 . . . R, 



A = ^2(r\e- pH A\r) / J^(r|e 



■I3H 



r . 



r=l 



r=l 



From equations ( |3.33| ) and (3.34) it follows that the statistical error is even reduced, 



A = (A) + C(l/v RZ). 



(3.35) 



(3.36) 
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For a general A, not commuting with H, we have to consider also the contribution of offdiagonal 
terms in the equation ( 3,32| ). Since the phases of random coefficients f3 rn are randomly distributed, 
we can expect that vectors (3 rn are approximately orthogonal, 



R 



R 



E ftrnPrm ~ J-j—(8 n m + Cnm/ V~R), 



(3.37) 



r=l 



where \Cnm\ = The relative contribution of offdiagonal terms is then given by 



£ (n|e-^ H A|m)C n 



n,m— 1 
n^vn 



5^(n|e-^A|n) 



(3.38) 



n=l 



It is not easy to estimate the ratio w in general. First we note that at (5 — ► we could choose |n) 
to diagonalize A, so offdiagonal terms could be avoided anyhow (for a static operator A). This is 
however not the case for low T. If we assume that the sign of £ nm is random and uncorrelated with 
matrix elements (n\e~^ H A\m), we can put an upper bound w(R, N s t) < \j\[R. There are however 
several arguments, e.g. operators A are usually local leading to a sparse-matrix representation, 
which seem to indicate a much smaller contribution of offdiagonal terms. 



To conclude, taking into account all assumptions mentioned, the approximation A ( 3.35 ) should 
therefore yield a good estimate of the thermodynamic average (A) at all T. For low T the error 
is expected to be of the order of 0(l/y/~R), while for high T the error is expected to scale even as 
0(l/\/N s tR). Since arguments leading to these estimates rely on several assumptions which are 
not easy to verify, it is essential to test the method for particular cases. 



3.6 Implementation and tests 

We comment now on the practical implementation of the FTLM and present few tests in order to 
get a quantitative estimate of approximations mentioned above. First we consider the calculation 
of static quantities. Joining the HTE and the random sampling we approximate the average of the 
operator A as 

at R M 
r=l j=0 

r=l j=0 

The sampling is over R random states \r) = \4> T ), which serve as initial functions for the M-step 
Lanczos procedure, resulting in M approximate eigenvalues with corresponding eigenvectors \4>j}- 
For a general operator A the calculation of eigenfunctions and corresponding matrix el- 
ements (ipj\A\r) is needed. On the other hand, the calculation effort is significantly reduced if 
[H, A] = and A can be diagonalized simultaneously. Then 

R M 
r=l j=Q 

In this case the evaluation of eigenfunctions is not necessary since the element (r\ipj) = Vj , equation 
(|3.4j ), is obtained directly from eigenvectors of the tridiagonal matrix H r M . 
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Few remarks on the implementation are in order here. Already in usual g.s. Lanczos calculations 
the use of symmetries of the model Hamiltonian is crucial in order to reduce computational and 
storage requirements. This is even more important when using the FTLM where the computational 
burden is increased due to the sampling and due to the calculation of matrix elements. At T > in 
general all symmetry sectors must be taken into account and these can differ significantly regarding 
the number of basis states they contain. Formulas (3.39), ( |3.40| ) must then be generalized to 



allow for the varying number of samples in each sector, so that sectors containing more states are 
more thoroughly sampled. If in the symmetry sector s containing Ng t basis states R s samples are 
evaluated, then the random sampling summation is modified as 

ifE-EfE- (3-) 

r=l s 4 r=l 

Usually we choose R s oc N* t . The number of Lanczos steps can also be taken as sector dependent, 
Ms < Ng t . This is important in sectors with small dimensions iVJ. 

Calculations on finite systems can be carried out on different lattices with various boundary 
conditions. For planar problems it is convenient to use tilted square lattices of sizes N = n 2 + m 2 
(Oitmaa and Betts 1978) with periodic boundary conditions (p.b.c). The translational invariance 
of lattice Hamiltonians is preserved on such systems, which makes crystal momenta k good quantum 
numbers and enables to reduce the basis of states and their dimension Ng t . 

Let us test the method for static quantities on the problem of a Heisenberg model on a two- 
leg ladder. We discuss the case J = J' (exchange equal along and perpendicular the ladder) 
which exhibits the spin gap. One of the most interesting quantities in this system is the uniform 
susceptibility xo = {(Si *) 2 ) /TN (defined and discussed in more detail for doped AFM in Sec. 6.2) 
and its T dependence. As a test we choose the model on a N = 2 x 8 ladder, for which the exact 



results obtained by full diagonalization are available (Barnes and Riera 1994). In Fig. 3.1a we study 
the influence of the number of Lanczos steps M on the accuracy of results. These are shown for 
fixed sampling R = 1028, while M is varied from 5 to 20. It is rather surprising that even with 
the smallest M = 5 in the largest symmetry sector we obtain very good agreement with the exact 
result not only at high T > J, but as well as at low T < J. In this case this is likely to be due to 
the gap in the energy spectrum, as expected for ladders with even number of legs. 



In Fig. 3J-b we fix M = 30, while the number of random samples R varies. Note that using 
the translational symmetry and the conservation of total S 1 " 1 = S z the maximum number of states 
in a symmetry sector is iVJ = 820, while the total number of states is N s t = 2 16 = 65536. The 
number of samples within each symmetry sector R s is chosen to be proportional to the number 
of basis states in the sector. In the largest sector Ng t = 820 this amounts to R s = 2 — 13 for 
the cases with total R = 180 — 1028, respectively. This corresponds to the sampling in the range 
R/Ngt = 0.003 — 0.016. We first observe that almost regardless of the sampling results agree closely 
with the exact one at higher T > J, as expected from the discussion of the statistical error in the 



equation ( 3.36; ). At T < J results start to disagree, still R s 3> 1 leads also to an improved accuracy 



in this regime. 



Let us turn now to the calculation of dynamical quantities. By joining equation (3.29) with the 



random sampling ( 3.35 ), we get the frequency-dependent correlation function 



R M M 



W)A) U « ^EEE^^ - e i + ^)m)mB\^)(^\A\r). (3.42) 



=1 j=o j=0 
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(a) 



(b) 



Figure 3.1: Uniform susceptibility xo vs - T for the 2 x 8-site Heisenberg spin ladder, as obtained 
using the FTLM with a) fixed sampling R = 1028 and various number of Lanczos steps M, and 
b) fixed M = 30 and different number of random samples R. Exact results are taken from Barnes 
and Riera (1994). 

The sampling is over R random states |r), resulting in M approximate eigenfunctions \ip\), |Vjf)> 
and corresponding e[,ej, respectively. 

At the full sampling R = N s t and within the chosen system the number of Lanczos steps M 
determines the number of exact frequency moments ///, analogous to g.s. moments ( gig ). This 
is evident from the expansion ( |3.28D at least for j3 — ► 0, while for lower T we are dealing with a 
double expansion, i.e. f3 and t series, and combined moments are exact. It follows that at least at 
high T the frequency resolution in spectra is Alu ~ AE/M, AE representing typically the energy 
span of the model. Since the information content in higher moments is limited, in particular due 
to finite-size effects, there is no point in using very large M, hence we restrict our calculation in 
most cases to M < 200. The effects of using a reduced sampling R <C N st are expected to be most 
pronounced at low T where moreover only a minority of symmetry sectors with the lowest energies 
contributes. 

Several tests for dynamical quantities within the t-J model (|2.6| ) have been already presented 
by Jaklic and Prelovsek (1994a, 1995c). Here we consider in addition the spectral function of a 
single hole injected in the undoped AFM, A(k,u) oc Re(ct (^^(O))^, defined and discussed in 

more detail in Sec. 7.2. We choose the system of 4 x 4 sites with J/t = 0.5 and k* = (tt/2,tt/2) 
corresponding to the g.s. wavevector. We compare T = results, being the most stringent test 
for the FTLM, with the g.s. ED results (Stephan and Horsch 1990, Eder et al. 1994). In Fig. |J 
we first show the convergence of the spectral function at T = with the sampling R s (note that 
for T = only g.s. symmetry sectors have to be considered) for a fixed number of Lanczos steps 
M = 180. Note that M determines also the number of correct T = frequency moments. We 
observe that the position of peaks in the spectrum is mainly unaffected by the sampling. Low- 
u> peaks are known to be quite accurate within the T = ED method, hence also their positions 
within the FTLM. Their intensities are less reliable at smaller sampling and also frequency moments 
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are expected to have larger errors. However, by increasing R s the accuracy of peak intensities is 
improved. 
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Figure 3.2: Spectral function A(k*,u) of a single hole in an AFM at T = 0. The g.s. ED result 
and the FTLM results for different random sampling R s are shown, at fixed M = 180. 



We investigate further the effect of reducing the number of Lanczos steps M. In Fig. 3J3 
the spectral function at T — » oo is calculated with R = 32 and varying M. We observe regular 
oscillations for lowest M = 30,60, appearing in frequency intervals Aw ~ AE/M, where AE is 
the maximum energy span in the model. We have only a partial explanation of this phenomenon, 
typical for high T. While the Lanczos algorithm obtains correct lowest and highest eigenvalues, 
Lanczos eigenvalues in the middle of the spectrum do not have any correspondence with true ones 
(as evident already from the discrepancy in their number M <C N^ t ) and appear almost equidistant. 
At T — ► oo they all contribute and yield observed oscillations. Since oscillations do not contain any 
relevant information, they can be easily smoothened out by a suitable filtering. They also become 
much less pronounced at lower T, where the predominant contribution is given by transitions from 
the states in the lower part of the spectrum to excited states. 



3.7 Finite size effects 

We introduce and justify the FTLM as a method to calculate T > properties on small systems. 
We also argue that choosing appropriate M, and the sampling R one can reproduce exact results 
to prescribed precision on a given system. In this sense the method is very effective, in its compu- 
tational effort comparable (although more time and memory consuming) to g.s. ED calculations 
on the same system. Still the well known deficiency of the ED method is the smallness of available 
lattices. Hence it is important to understand the finite size effects and their role at T > 0. In the 
following we predominantly study planar systems corresponding to the tilted square lattice with 
p.b.c. (Oitmaa and Betts 1978) where N = n 2 + m 2 . Mostly, we employ N = 16,18,20,26, as 
presented in Fig. |3.4| . 
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Figure 3.4: Tilted square lattices with p.b.c. corresponding to different sizes. 



We claim that generally T > reduces finite size effects. This is related to the fact that at 
T = both static and dynamical quantities are calculated only from one wavefunction |\Po)> which 
can be quite dependent on the size and on the shape of the system. In particular, g.s. spectra of 
dynamical quantities (see Dagotto 1994), e.g. the optical conductivity (Sega and Prelovsek 1990) 
and the single-particle spectral function (Stephan and Horsch 1991), quite generally appear as a 
restricted number of delta functions. While lowest frequency moments, in the sense of equations 
( 3.1S| ), can be quite representative of a large system, the peak- like structure and details of spectra 
are mostly not. 

T > introduces the thermodynamic averaging over a larger number of eigenstates. This 
reduces directly finite-size effects for static quantities, whereas for dynamical quantities spectra 
become denser. From the equation (3.42) it follows that we get in spectra at elevated T > 
generally RM 2 different peaks leading to nearly continuous spectra. This is also evident from 



high-T result in Fig. |3.3|, as compared to the T = result in Fig. 3.2 



The effect of T > can be expressed also in another way. There are several characteristic length 
scales in the system of correlated electrons, e.g. the AFM correlation length £, the transport mean 



21 




-20 



(a) 



(b) 



Figure 3.5: Many-body levels for: (a) the 2D Heisenberg model on N = 16 sites, (b) the 2D t-J 
model with = 2 on iV = 10 sites. In both cases only the S z = 0, k = sector is presented. 



free path l s , etc. These lengths decrease with increasing T and results for related quantities have a 
macroscopic relevance provided that the lengths become shorter than the system size, e.g. l s < L 
where L is the linear size of the system. This happens for particular T > T s , where clearly T s 
depends also on the quantity considered. For certain quantities one can monitor such conditions 
directly, e.g. for l s as discussed in more detail in Sec. 5. 

As a criterion for finite size effects we use the characteristic finite-size temperature TV S . It is 
chosen so that in a given system the thermodynamic sum 

Z(T) = Tre~^ H ~ Eo) (3.43) 

is appreciable, i.e. 2{Tf s ) = Z* 3> 1. 

To get the size dependence of Tf s it is important to understand general features of many body 



spectra. In Fig. 3.5a we present the levels of the Heisenberg model on a square lattice with N = 16 



sites for the S z = 0, k = sector. Note that energy span is AE oc N J while the number of 

This means that the density of states far from spectral edges scales 



•>N 



"st 



exponentially with N. On the other hand, near the edges spectra become sparse and the spacing 
between lowest levels decreases rather slowly with the size, i.e. Ae oc N~ p . We expect that also 
Tj s scales as Ae. Still the character and the density of low-lying states can change qualitatively 
from one regime to another. In Fig. [3.5| b we show for comparison the levels within the t-J model 
with Nh = 2 holes on the tilted square lattice with N = 10 sites, again only for the S z = 0, k = 
sector. While both cases have similar N^., it is evident that the low-energy regime shows much 



higher density of states in doped t-J model, Fig. |3.5[ b, indicating a large degeneracy of states at 
The FTLM is best suited just for quantum many-body systems with a large degeneracy of 
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Figure 3.6: The variation of Tf s with doping Ch in the t-J on N = 18 sites with J/t = 0.3. Curves 
correspond to different choices of Z* = Z(Tf s ). 

states, i.e. large Z at low T. This is the case with doped AFM and the t-J model in the strong 
correlation regime J < t. To be concrete we present in Fig. |3.6| the variation of Tf s with the doping 
Ch = Nh/N, as calculated from the system of N = 18 sites and J/t = 0.3. For convenience we fix 
Tf s with criteria Z* = 30, 100, 300, respectively. It is indicative that Tf s reaches the minimum for 
intermediate (optimum) doping Ch = c* h ~ 0.15, where we are able to reach Tf s /J < 0.3. Away 
from such optimum case Tf s is larger. In the undoped (and underdoped) AFM this happens due to 
rather large finite-size gaps in magnon excitations, while in the overdoped system electrons behave 
closer to free electrons with a nearly unrenormalized bandwidth and hence large gaps between 
single-electron excited states. We claim that small Tf s and related large degeneracy of low-lying 
states are the essential features of strongly correlated system in their most challenging regime, being 
a sign of a novel quantum frustration. On the other hand this gives an advantage to the FTLM 
which performs best where several other methods, like QMC, fail due to the same frustration (sign) 
problems. 

3.8 Relation to other numerical methods 

It is not our intention to give an exhaustive overview of other numerical (or partly analytical) 
methods which are used in the analysis of the problem of strongly correlated electrons, and for 
doped AFM in particular. We mainly list below methods which are alternative to the FTLM, 
together with their limitations and advantages. Since T > calculations are often used as a proper 
approach to T = results, as is also the case for several quantities evaluated within the FTLM 
furtheron, we comment also on some g.s. calculations. 

Clearly the closest relation is with the ED studies of small systems. Apart from the FTLM only 
few studies of T > properties of the t-J model and of the Heisenberg model have been performed 
(Tohyama et al. 1993, Sokol et al. 1993, Tsunetsugu and Imada 1997), restricted to small ./V and in 
particular small Ng t < 1000 due to the full diagonalization employed in these calculations. We note 
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that such restricted N s t can influence in particular dynamical spectra, which can show up finite 
size effects even at higher T. We should also note that quite an analogous method to the FTLM, 
based instead on the Chebyshev iteration, has been introduced recently (Silver and Roder 1994), 
but has not been exploited much so far. 

Much more extensive ED calculations have been performed for g.s. properties (see Dagotto 
1994). While systems considered could be possibly somewhat larger than those reachable within 
the FTLM, we claim that in particular dynamical spectra as evaluated within the g.s. should be 
interpreted with care due to their sparse structure. This seems to be the case with the g.s. optical 
conductivity, the spin structure factor, spectral functions etc., as discussed in detail furtheron. 

Closely related to the FTLM (see Sec. 3.3) is the HTE approach, which in principle deals with 
a large (infinite) system and is one of few straightforward (at least partly) analytical methods for 
correlated systems. It uses (3 as a small parameter for the series expansion, while appropriate 
extrapolations (e.g. using Pade approximants) are needed to obtain results for T in the physically 
interesting regime. So far it has been used in several studies of the t-J model, e.g. to address 
the question of the phase separation (Putikka et al. 1992), magnetic correlations (Singh and 
Glenister 1992a), the momentum-distribution function (Singh and Glenister 1992b), the charge- 
spin separation (Putikka et al. 1994), the Hall effect (Shastry et al. 1993) etc. The advantage of 
the method is the absence of any finite cluster bounds. On the other hand the method requires a 
careful and a nonunique extrapolation procedure which is reliable only for static quantities. 

Widely used QMC technique yields in general results for T > 0. There are various methods, 
which have been covered in several recent reviews (von der Linden 1992, Suzuki 1993). We mention 
here only approaches which are relevant for studies of planar undoped and doped AFM. The world- 
line QMC has been very successful in the evaluation of static properties of the Heisenberg model 
(see Manousakis 1991). Analogous results have been obtained via the QMC for the insulating 
Hubbard model at half-filling (Hirsch 1985). Away from the half-filling the sign problem becomes 
the major difficulty for the QMC studies of fermionic models. It is particularly severe at low T 
within the intermediate-doping regime, being thus connected with the large degeneracy of fermionic 
states. Still various static quantities have been evaluated within the Hubbard model as a function 
od doping both for T = and T > (see Dagotto 1994). 

The calculation of dynamical quantities within the QMC is possible via the deconvolution of the 
imaginary-time dynamics into a real-frequency one using the maximum entropy analysis (Jarrell et 
al. 1991). The latter appears to be quite delicate due a large influence of statistical errors. Still 
there has been in recent years several studies of dynamic properties of spin systems (Makivic and 
Jarrell 1992) and of the planar Hubbard model, in particular of spectral functions (Bulut et al. 
1994, Preuss et al. 1995, 1996). In spite of much larger systems reachable within QMC studies, it 
is well conceivable that due to inherent difficulties QMC results for dynamics are less reliable than 
those obtained within the FTLM. 

There are other powerful numerical methods which are only partly relevant to studies of doped 
AFM. Particularly successful and promising is the Density Matrix Renormalization Group (DMRG) 
approach, as developed by White (1992) and extensively applied to problems of correlated electrons. 
While designed mainly for ID systems, it has been extended to ladder systems (White and Scalapino 
1997a) as well as to planar models (White and Scalapino 1997b). So far appropriate generalizations 
in order to study the T > and dynamical properties have only been attempted (Pang et al. 1996). 
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4 Thermodynamic properties 



We first consider thermodynamic properties of the t-J model (Jaklic and Prelovsek 1996). These 
include quantities directly derivable from the grand-canonical sum f2: free energy density J-, chem- 
ical potential /J,, charge compressibility k, entropy density s, specific heat Cy etc. Some of these 
have been already studied using other methods. Within the HTE some thermodynamical quantities 
have been calculated within the t-J model. Results indicate the ferromagnetic phase at J <C t and 
at low doping (Putikka et al. 1992), but also a large enhancement of the entropy in a doped AFM 
(Putikka, unpublished). Within the Hubbard model the projector (T = 0) QMC method has been 
employed to study the phase diagram and to calculate the charge compressibility (Furukawa and 
Imada 1992). Several calculations of the chemical potential vs. doping within the Green's function 
QMC have been presented in order to establish the regime of the phase separation (Kohno 1997, 
Hellberg and Manousakis 1997), although with contradictory conclusions for the most interesting 
regime J < t. 

In order to study continuously varying particle densities, we perform the averaging within the 
grand-canonical ensemble, involving all possible numbers of electrons N e , 

(A) = Y,^N e [e-^ H -^A]/Y,^N £ e-^ H -^\ (4.1) 

where [i is the chemical potential. For each N e the problem thus reduces to the evaluation of the 
canonical thermal average, which we achieve with the FTLM as described in Sec. 3. 

The implementation of the FTLM can be further simplified for operators A which are conserved 
quantities, i.e. commute with H, as shown in Sec. 3.6. Examples include H itself, the particle 
number N e , the total spin S z etc. By choosing random functions |r) to have good quantum numbers 
A e and S z , we can evaluate the expectation value of an arbitrary function f(N e , S z , H) 

(/) - ^EEKri^)i 2 /(^,5r, e pe-^-^), 

r=l j=0 

n - ^EEKrl^lV^-^. (4.2) 

r=l j=0 



As noted in the equation ( 3.4C| ), in this case |(^|'i/'j)| 2 can be evaluated from the tridiagonal matrix 



directly. Since also M < 100 is enough, the reorthogonalization of Lanczos functions can be avoided. 
This eliminates the need to store wavefunctions |^>p and systems with considerably larger N s t can 
be studied. Consequently the computational effort in this case is equal to that of a g.s. Lanczos 
procedure, repeated R times. We employ in the following typically R ~ 200 — 1000 in each N e 
sector. Calculations are performed on systems with N = 16, 18, and 20 sites for arbitrary filling 
c/i, while for the undoped system we reach A^ = 26 sites. Note also that we fix J/t = 0.3. 

4.1 Chemical potential 

We first analyze the hole chemical potential /x^, = — fi as a function of T and of the hole density 



Ch = 1 — (N e )/N. Results are obtained by first calculating at fixed /i,T from equations (4.2) 



with / = N e , and then inverting the dependence Ch(^,T). In Fig. [LlJ we present curves Hh(T) for 



several c%. Note that for thermodynamic quantities results seem somewhat less sensitive to finite- 



size effects, hence we follow them in Fig. |4.l| to T ~ 0.05 t (for Ch < 0.1 we would still estimate 
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Figure 4.1: Hole chemical potential Hh vs. T at several dopings c^. 



higher Tf s ~ 0.1 t). In order to interpret Hh(T) at low <C 1 it is essential to note that at T = 
the system contains no holes in the equilibrium provided \Xh < /J^. For chosen J/t = 0.3 it has been 
established ~ —1.99 t (Dagotto 1994), related to the minimum energy of a single hole added to 
the undoped AFM. 

Analyzing Fig. 4.1 we mostly do not find a T 2 dependence of fx^ at low T, as expected for a 



normal LFL, except within the extremely overdoped regime Ch > 0.3. In particular, in a broad 
range 0.05 < < 0.3 we find a very unusual roughly linear variation, 



Hh(T) = n h (T = 0) + aT, T fs <T< J, 



(4.3) 



whereby the slope a changes the sign at c/j = c* h > 0.15. It is remarkable that the marginal doping 
c h appears to be quite system independent, as checked quantitatively for different system sizes 
TV = 16 — 20. 

The marginal c* h shows up again in Fig. |4.2| , displaying the variation c/j(A/i) at various temper- 
atures T. Here A/j, = \i + ^ h is the difference to the undoped AFM case and is displayed in eV to 
allow the comparison with experiments, using the usual correspondence t = 0.4 eV. c* h represents 
in this case the crossing of curves at different T, whereby A//(T) is essentially pinned at the value 
Ajj, ~ —0.17 t. This pinning is active in a wide range of T. Analyzing the regime < (? h in 
Fig. 4.2, we note again that c^(T — ► 0) remains finite only for fj,^ < /x^. One can evaluate from 
these results the compressibility of the hole fluid k oc —dch/dfj,. We find that k < oo for all T > 0, 
indicating the absence of the phase separation in the system at chosen J/t. This is in contrast with 
some recent QMC studies (Hellberg and Manousakis 1997) claiming the phase separation in the 
t-J model at all J/t, as first put forward by Emery et al. (1990). Our results in Fig. ^ do not 
support such a behaviour, at least not in the range T/t > 0.05. Still Fig. 4.2 reveals that n{T — > 0) 
is increasing and becoming larger on approaching A/i ~ 0. 

A proper interpretation of the regime Ch — > at low T is clearly one of the major challenges. 
One frequently used picture is that holes doped in an AFM could be described as degenerate 
fermions with small FS - hole pockets (Trugman 1990, Eder and Becker 1991). In this case one 
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Figure 4.2: Hole concentration vs. Afi at various T. For comparison we show experimental 
results for LSCO at low T, obtained from the shift of photoemission spectra by Ino et al. (1997a). 



would expect in a 2D lattice k = N(n) = l/2nt*, where M(e) is the single-particle (hole) density 
of states (DOS) and t* an effective hopping parameter. Since the hole effective mass can be quite 
enhanced, i.e. f « t, k can become quite large. Results in Fig. O put a lower bound to the 
possible enhancement, i.e. t/t* > 5. 

Analyzing T = QMC results for the Hubbard model near half filling a variation Ch(^h), 
similar to ours, has been found (Furukawa and Imada 1992, Assaad and Imada 1996). Results were 
interpreted in terms of a singular behaviour k oc (fih — A^) -1 ^ 2 - I n contrast to the hole-pocket 
picture, the latter form does not allow for a regime with a degenerate gas of holes, at least not 
with holes having a doping-independent mass. From our results it is hard to exclude any of these 
scenarios, whereby the regime of hole pockets should be in any case restricted to very low doping 
c h < 0.1. 

It is tempting to interpret the existence of the marginal concentration c* h as a change of the 
character of the FS. To establish the relation, we have to rely on arguments which apply to the 
gas of noninter acting fermions. A simple Sommerfeld expansion yields that the electron density at 
fixed (jl is given by 

c e (T) = c e (T = 0) + (7rfe * r) V (M). (4.4) 

6 

Indirectly this gives an information on the FS, since one would plausibly associate A/"'(/i) > for 
c e < 1 with a large electron FS, and oppositely M'{n) < with a hole- like FS or small hole pockets 
vanishing for c e — > 1. At least for free electrons it is easy to establish the connection of M'(fi) with 
the curvature K~ l of the FS (Jaklic and Prelovsek 1996), analogous to the relation for the Hall 
resistivity (Tsuji 1958). At least in the region of the k space, where the effective-mass tensor is 
positive-definite, M'{p) > implies also that the average FS curvature K~ x is positive. 

The observed nonquadratic T dependence in Fig. [4.1| questions the interpretation in terms of the 



free-electron DOS ( |4.4|) . Still we may interpret dch/dT < for Ch > c* h , as deduced from Fig. 4.2 , 
as an indication for Af'([i) > 0, i.e. positive average curvature of the FS. This in turn implies a 
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(a) (b) 
Figure 4.3: Entropy density s vs. T for different dopings Ch for a) J/t = 0.3, and b) J = 0. 

transition at Ch ~ c£ from a hole-pocket picture at low doping (Trugman 1990, Eder and Becker 
1991), to an electron-like large FS (Stephan and Horsch 1991, Singh and Glenister 1992b). 

Recently the variation of [i with the hole doping in LSCO has been deduced experimentally 
from the shift of photoemission spectra by Ino et al. (1997a). For comparison we plot also these 
results in Fig. [O^, noting that they apply to low T in terms of our model parameters. The overall 
agreement is quite reasonable, taking into account the uncertainty of PES results. Again the flatness 
of n(ch < c* h ) is quite remarkable, indicating on the possibility of divergent k — > oo for Ch — > 0. As 
noted already by the authors the variation /i(c^) is highly nontrivial and cannot be accounted for 
by a simple LFL results as e.g. obtained in band structure calculations. 

4.2 Entropy 

Let us consider the entropy density (per unit cell) 

,4( teln!i+ M^W), (4 . 5) 

where averages and fi are calculated using the FTLM (Jaklic and Prelovsek 1995b, 1996) with 
equations (|4.2| ). Within the t-J model s has been studied also via the HTE (Putikka, unpublished), 
while the QMC method has been recently used to calculate the entropy within the Hubbard model 
(Duffy and Moreo 1997). 

The T variation of s at various is presented in Fig. [4.3ja. Note again that within the grand- 
canonical calculation Ch can be followed continuously. Results shown for N = 20 are quite close to 
those for lattices with N = 16, 18 sites provided that T > 0.1 t. It is evident that in the undoped 
AFM s(T) at low T < J/2 is consistent with the magnon contribution s oc T 2 . This dependence 
changes however already for smallest finite doping = 0.05 to s oc T a with a < 1. 

In order to understand the role of AFM correlations induced by J > on the entropy s, we 
present in Fig. [4.3j b also results obtained for J = 0. It is clear that here the undoped case = 
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Figure 4.4: s vs. Ch at several T, calculated in a system with N = 20 sites. We show for comparison 
also experimental results for LSCO by Loram et al. (1996) at highest T = 320 K ~ 0.07 t. 



is singular, since s(T > 0) = In 2. It is also plausible that for Ch > 0.3 we find s(T) essentially 
independent of J, i.e. the spin exchange becomes irrelevant in the overdoped regime. Comparison 
again confirms that the role of J is crucial in the underdoped regime and at optimum doping. 



Alternatively we can discuss the doping dependence of s, as shown in Fig. |4.4| at different T < J. 
As realized already from the discussion of the thermodynamic sum Z(T) and related Tf s (ch) in 
Fig. |3.6| , the entropy displays a broad maximum at c\ t ~ 0.15, indicating the highest density of 
many-body states in the optimum-doping regime. The appearance of the maximum in s(c/J is 
intimately related to fi(T) discussed in Sec. 4.1. Namely from general thermodynamic relations 
(equality of mixed derivatives) for the free energy density T{ch,T) it follows 



ds 
dc h 



dT 



taking into account that s = —dJ-jdT and Hh = dJ- /dch- The relation (4.6) connects s(ch) 



(4.6) 



max 



with the pinning of Hh{T) seen in Figs. 4.1 , [4.2| at the optimum doping Ch ~ c* h . 

Besides the enhancement of s with doping, the surprising fact is also its magnitude at T < J, 
i.e. at small T in terms of model parameters, s in the optimum regime appears very large, e.g. at 
T = 0.1 t = J/3 the entropy per site is s ~ 0.39 ks, which is almost 40% of s(T = oo) for the same 
Ch, although T < J and the energy span of excitations extends well beyond the scale AE > t. This 
should be contrasted with the situation in an undoped AFM, where s becomes relatively significant 
only for T > J/2, and saturates for T > J. Moreover, in the case of noninteracting fermions one 
gets s ~ ks only at the Fermi temperature T% ~ W, where the bandwidth is W = 8 t. On the other 
hand, by introducing the degeneracy temperature within the t-J model as s(r^ e3 ) = s(T = oo)/2, 
we get for ~ c* h only T^eg ~ 0.17 t, being small in comparison with any reasonable effective QP 
bandwidth. 

It is indicative that the entropy of such a magnitude has been deduced from the electronic 
specific heat measurements in oxygen deficient YP^CusC^-^ (YBCO) materials (Loram et al. 
1993). E.g., for the optimally doped material with 5 = 0.03 at T = 300 K the experimental 
result is As = 0.35 ks per planar copper site (As = 0.7 ks per formula unit), relative to the 
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Figure 4.5: Specific heat Cy vs. T for the undoped AFM for several system sizes. 



undoped 5 = 1 sample. We find the corresponding value As = s(ch = 0.15) — s(c^ = 0) ~ 0.30 ks 
at T = 0.1 t ~ 450 i^. Recently, s has been measured also for LSCO in a large doping range 
< x < 0.4 (Loram et al. 1996). We plot results at fixed T = 320-fT ~ 0.07 t as a function of doping 
in Fig. |4.4j , for comparison with our model results. The qualitative and quantitative agreement is 
quite promising, also in view of possible uncertainties in the experimental determination of s. When 
comparing results we note that our curves start at higher T and so the main difference is in the 
location of the entropy maximum, which in LSCO appears at somewhat higher Ch ~ 0.22. 



4.3 Specific heat 

The same results can be discussed in terms of the specific heat (per unit cell) 



(4.7) 



which can be as well represented directly with expectation values at given T, analogous to the 
equation ( fLq) and a differentiation with respect to T is not needed. Let us first show as a test 
Cy for the undoped AFM (Heisenberg model) for several system sizes. In Fig. |4.5| results are 
shown for systems ranging in size from 16 to 26 sites. Except at the lowest T < J/3, results 
do not vary appreciably with the system size, particularly regarding the position and the height 
of the maximum. Our results seem to be even superior to those obtained by the QMC method 
(Gomez-Santoz et al. 1989). Calculated Cy is strongly T dependent in the whole T range, with a 
maximum at T ~ 2 J/3, and as expected Cy oc T 2 at low T consistent with the magnon excitations 
dominating this regime (Manousakis 1991). 

In Fig. |4.6| we present Cy(T) at different c^. As the AFM is doped, Cy(T) still exhibits a 
maximum, which is however strongly suppressed and gradually moves to lower T with increasing 
Ch- The peak can be attributed to the thermal activation of spin degrees of freedom. The latter 
are still characterized by the exchange scale J which persists in the doped system, as observed 
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Figure 4.6: Cy vs. T for different hole dopings c^. 



also in dynamical spin correlations (Jaklic and Prelovsek 1995a) discussed in Sec. 6. The exchange 
energy scale however disappears in the overdoped regime Ch > 0.3. Results indicate a possible LFL 
behaviour with Cy ~ 7T only for T < 0.1 t. It is characteristic (and consistent with the vanishing 
role of J) that in the optimally doped regime ~ 0.2 we find Cy{T) ~ const for 0.15 < T/t < 1, 
being far from a FL behaviour. 

Results in Fig. [4.6| confirm the recent conjecture (Vollhardt 1997) that in correlated systems 
the specific heat Cy(T,X) can show universal crossings as a function of a thermodynamic variable 
X. In our case we consider X = Ch, and realize that Cy(T) cross for different Ch at two T, whereby 
the lower crossing at T ~ 0.13 t seems to be nearly independent of c^. 

5 Electrical Properties 

The anomalous normal-state character of electrical transport properties has been realized since the 
discovery of high-T c cuprates and remains the challenge for theoreticians ever since. Among these 
properties the prominent example is nearly linear in-plane resistivity p oc T in the normal state 
(for a review see lye 1992, Batlogg et al. 1994). It is however an experimental fact that such 
a behaviour is restricted to the optimum doping regime, while deviations from linearity appear 
both in the underdoped and overdoped regimes, being still universal for a number of materials 
with a similar doping (Takagi et al. 1992). The d.c. resistivity is intimately related to the optical 
conductivity (t(u>), which has been also extensively studied (for a review see Tanner and Timusk 
1992) and shows in the normal state the unusual non-Drude behaviour (Schlesinger et al. 1990, 
Romero et al. 1992, Cooper et al. 1993, El Azrak et al. 1994, Puchkov et al. 1996, Startseva 
et al. 1997). Another challenging set of experimental findings concern the d.c. resistivity p c {T) 
perpendicular to the Cu02 planes and the corresponding a c (u) (see Uchida 1997), which we will 
not consider here. 

The central question is whether these anomalous static and dynamical transport properties can 
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be accounted for by strong correlations alone, or possible other mechanisms such as the electron- 
phonon coupling have to be invoked (Zeyher 1991). In spite of considerable efforts so far there 
are very few microscopical theories of electron transport dealing with planar or higher-dimensional 
strongly correlated systems. 

Brinkman and Rice (1970) solved the problem of a single mobile hole in the extreme case 
J = within the retraceable-path approximation (RPA) and evaluated the d.c. mobility hq(T). 
Analogous results have been obtained via the HTE by Ohata and Kubo (1970). Within the RPA 
also a(uj) has been evaluated (Rice and Zhang 1989). It is not easy to find the range of validity and 
relevance of these results, in particular for J > where one would expect possibly a crossover to a 
different relaxation mechanism for T < J. Nevertheless it is clear that the analysis, treating holes as 
independent, is more appropriate for the weak-doping regime Ch <C 1, at least when the behaviour 
at low uj < t and T < t is concerned. It appears much more difficult to approach analytically the 
electron transport at Ch > 0. An attractive proposal remains that of spinons and holons as basic 
low-energy excitations (Anderson and Zou 1988), as well as related gauge theories (Nagaosa and 
Lee 1990) and slave-boson approaches, which have been applied also to the calculation of optical 
conductivity (Bang and Kotliar 1993). Very fruitful have been recent studies of infinite-dimensional 
models, in particular for the Hubbard model (for a review see Pruschke et al. 1995, Georges et al. 
1996), which allow also for the evaluation of ct(uj) and p(T). Since these numerical results are also 
hard to interpret, it is still under debate to what extent they contain features relevant to lower 
dimensions, e.g. to planar systems discussed here. 

Several conclusions on the charge dynamics have been reached using the ED dealing with the 
g.s. behaviour. For a single hole in the AFM the planar optical conductivity ct(uj) (Sega and 
Prelovsek 1990, Poilblanc et al. 1993) and the charge stiffness (Zotos et al. 1990) have been 
interpreted in terms of partially coherent hole motion with a substantially enhanced effective mass 
(at J/t < 1) and with a mid-infrared peak at uj ~ 2 J. Analogous results were presented for larger 
doping (Dagotto 1994), however with considerable finite-size effects, so that their interpretation 
does not appear well settled. We discuss in the following results for the charge transport in the 
t — J model as obtained with the FTLM, in part already presented elsewhere ( Jaklic and Prelovsek 
1994b, 1995c). 

5.1 Current response 

Let us consider the real optical conductivity cr(uj), which is a tensor in general. We are dealing here 
with a system without any external magnetic field. On a square lattice the tensor is then diagonal, 
so that a aa (uj) = a(uo). Within the linear-response theory (see e.g. Mahan 1990) the regular part 
of a(uj > 0) is given by 

1 — p-l 3 ^ 1 r°° 

a reg (uj) = el C(uj), C (uj) = -Re dte^{j a (t)j a (0)), (5-1) 

uj i\ j 

where j is the (total) particle current operator. In a finite system one can write C{uj) in terms of 
exact eigenstates |* n ) with corresponding energies E n , 

C ^ = J^T, e-^\(^ m \j a \^ n )\ 2 6(co -E m + E n ). (5.2) 

It is however well known that in general one has to take into account also the singular contribution 
to the charge dynamical response, i.e. 

a(uj) = 2irelD c 5(uj) + a reg (uj), (5.3) 
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where D c represents the charge stiffness. We study in the foliowing cr(u) in more general tight- 
binding models, e.g. including also the n.n.n. hopping. Since the analysis of this case, together 
with the derivation of a proper D c and the optical sum rule, is not usual in the literature we present 
it shortly below. 

We follow the approach by Kohn (1964) introducing a (fictitious) flux cf) through a torus repre- 
senting the square lattice with p.b.c. Such a flux induces a vector potential A, being equal on all 
lattice sites. In lattice models with a discrete basis for electron wavefunctions A can be introduced 
with a gauge transformation (Peierls construction) cj s — ► cj s exp(— ieoA ■ Rj), which effectively 

modifies hopping matrix elements. Taking A as small we can express the modified tight-binding 
Hamiltonian allowing also for more general hopping elements 



. Lint 



« H(0)+e A-j+-V-A-rA, (5.4) 
where R\j = Rj — Ri, r is the kinetic stress tensor, and 

3 — i ^ ~] tjjRjjCjgCjsi 
i,j,s 

t = ^^tijRij (£> RijCj g Ci s . (5-5) 



Note that in usual n.n. tight-binding models r is directly related to the kinetic energy operator, 

The electrical current j e is from the equation fl5.4|) expressed as a sum of the particle-current 
and the diamagnetic contribution, 

J e = -8H/dA = -e j - elrA. (5.6) 

The above analysis applies also to an oscillating A(t) = A(uj)exp(—iuj + t). This induces an electric 
field in the system E(t) = —dA{t)/dt. We are interested in the response of {j e )(oj). Evaluating (j) 
within the standard linear response (Mahan 1990), and with A(u) = E(lo)/iuj + , we arrive at the 
complex optical conductivity 

= ^«T> -*(*)), 
poo 

XH = i / dte^ +t ([j(t)J(0)}). (5.7) 
J o 

Complex = ct{lo) + i&"(oS) satisfies the Kramers-Kronig relation. Since %(cj — > oo) — * 0, we 
get from the equation Q5.7|) a condition for o-"{uj — * oo), 



f 

J — ( 



2 



^ = ^W, (5.8 



which corresponds to the optical sum rule. It reduces to the well known one for continuum electronic 
systems, as well as for n.n. hopping models where (T aa ) = —(Hki n )/d (Maldague 1977). We can 
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now make contact with the definition (|5.3|). From the expression (|5.7|) it follows 
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^ = ^2 lim n w C(w) = 771T?[(^a) -Xaa(°)]" ( 5 - 



5.2 Charge stiffness 

Nonzero charge stiffness = D C (T = 0) > is a characteristic signature of a metallic state (Kohn 
1964, Scalapino et al. 1993), in contrast to an insulator with D® = 0. The evaluation of D® has 
been recently applied to a number of correlated fermionic systems, both analytically for ID systems 
(Shastry and Sutherland 1990) and numerically for planar Hubbard and t-J models (see Dagotto 
1994). At T > one would expect for normal resistors D c = 0. Since we are working with small 
systems with p.b.c, where ballistic response of carriers can persist at T > 0, we find in general 
D C (T) ^ 0. Note however that recently a nontrivial possibility of a nonergodic behaviour in a 
macroscopic limit has been realized, i.e. with D C (T > 0) > 0, being related to the integrability of 
the fermionic model (Castella et al. 1995, Zotos and Prelovsek 1996). 

The t-J model at half-filling has D C (T) = at any T, since charge fluctuations are projected out 
by construction of the model Q2,6j ). For a doped system one expects D® oc Ch at low Ch- Studying 
the charge transport at T > within the planar t — J model we adopt the view that there exist 
scattering processes which in a macroscopically large systems cause D C (T > 0) = 0, i.e. the model 
is ergodic. However, in our numerical calculations we are dealing with small systems and a small 
portion of the total current can propagate through the system unscattered, thereby establishing 
in the system a persistent current. Hence we characterize observed D C (T > 0) 7^ as a finite-size 
artifact. Nevertheless, a variation of D C (T) brings a valuable information. At T > the scattering 
processes result in a finite mean- free path of charge carriers l s (T)- When the linear system size L 
exceeds l Sl it is reasonable to expect that D C (T) ~ 0. On the other hand for L < l s (T) we obtain 
D C (T) ~ D®. By following D C (T) we can thus independently monitor the transport mean free 
path. Since l s is strongly T-dependent, we can approximately locate the crossover temperature as 
L ~ l s (T fs ). 

We can express D c in the t-J model on a square lattice from the equation (|5.9|) as 



1 ,„ , 1 



00 



D c = - — (H km ) 2 / a{cv)du. (5.10) 

Within the FTLM {H^in) and o~(oj) are calculated separately. It should be however mentioned that 
due to finite M there could be some ambiguity in the cutoff employed at low u, the problem which 
is however restricted to otherwise unproblematic T 3> t regime. 

Results for a single hole in the t—(J = 0) model, as obtained by Jaklic and Prelovsek (1995c) 
show that D C (T) interpolates quite smoothly between D c = at T = 00 and D®, the crossover 
becoming sharper in larger systems. For the AFM case results for D C (T) are less regular, as 



presented in Fig. 5.1 for various > 1 on a system with N = 16 sites. When judging the extent 
of deviations of D c from zero at T > 0, it is useful to compare values with the maximum possible 
ones, i.e. with the value of the sum rule at T = 0, D max = \{Hki n )(T = 0)|/4iV, as follows from 



the equation (|5.10| ). We notice from Fig. 5^ that D c typically shows a rather abrupt transition 



from D c ~ to D c 7^ at the crossover temperature Tf s , depending mainly on Ch- For T < Tf s 
the variation D C (T) can become quite unphysical, i.e. we get in some cases even D c < 0. These 
phenomena are influenced by particular p.b.c. and more sensible results can be obtained by the 
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Figure 5.1: Normalized charge stiffness D c /D max for various number on holes Nh in a system with 
N = 16 sites. 



introduction of twisted boundary conditions or fluxes (Poilblanc 1991). Nevertheless we are here 
interested only in the regime T > Tf s . It follows again that Tf s is minimum, i.e. Tj s ~ 0.1 t, for 
the intermediate doping 0.1 < Ch < 0.3. We should stress the striking message that at intermediate 
doping even at such low T, representing for cuprates T ~ 450K, the mean free path does not exceed 
l s ~ 4 sites and this entirely due to correlation effects. 



5.3 Single-hole mobility 

It is expected that at low doping the conductivity scales linearly with doping, hence it is meaningful 
to introduce the dynamical mobility /j,(u)) which is a single-hole property 

o(uj) = e c h fi(uj), (5.11) 

but still highly nontrivial for a Mott or an AFM insulator. 

The most conclusive theoretical results (in 2D or higher D systems) have been so far obtained 
for a problem of a single mobile hole introduced into a reference insulator. Brinkman and Rice 
(1970) solved the problem for J = within RPA. They pointed out of an essentially incoherent hole 
motion and evaluated the d.c. mobility fio(T), exhibiting hq ctT loiT > t. fi(uj) within RPA (Rice 
and Zhang 1989) shows an incoherent motion, resulting in a slow non-Drude fall-off fi(uj) oc 1/uj for 
larger uj. The RPA has been recently justified and applied more rigorously for infinite-D lattices 
(Metzner et al. 1992). An analogous approach is the evaluation of frequency moments of /u(w), 
starting at T » t, as applied to the J = problem by Ohata and Kubo (1970). On the other 
hand, fi(uj) at T = has been in recent years well established by numerical studies of small systems 
via the ED method (Sega and Prelovsek 1990, Dagotto 1994). Nevertheless there are important 
unsolved questions even for the single-hole problem. Is ^i(w) on a planar lattice qualitatively and 
quantitatively well described within the RPA, at least for T > t ? Which are new qualitative 
dynamical features at T < t, both for the J = and the J > case ? 
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Results for (jl(u) at J = 0, obtained via the FTLM by Jaklic and Prelovsek (1995c), show an 
overall agreement with the RPA. However, in contrast to the smooth RPA curve the actual (i{ui), 
evaluated at T > f, seems to be nonanalytical, i.e. it shows a cusp at u = 0. The phenomenon 
seems to be characteristic for J = 0, but not for J > 0. (l(uj) for J = retains its high-T form 
for all T > t, but changes qualitatively for T < t. Here the central peak due to the formation of 
the FM polaron (Nagaoka 1966) starts to emerge at low uj, and the incoherent broad background 
steadily vanishes on approaching T —* 0. 

We are more interested in the AFM case, as shown in Fig. For T > t the spin background 
is disordered and the mobility retains the high-T form (l{ui) oc (3, leading to 1/T. As seen 

from Fig. |5.2| , there is a qualitative change already in the regime J < T < t. Namely, it appears 
that here (i{ui) has a weaker T-dependence. For T —* we are however approaching a nontrivial 
response of an AFM polaron, which has been analyzed numerically by several authors. At T = 
one expects in (i{u)) for w->0a coherent response of an AFM polaron with an enhanced mass, but 
also a nonvanishing incoherent part at uj > J (Sega and Prelovsek 1990, Dagotto 1994). The latter 
component seems to have a nontrivial internal structure being related to the mid-infrared peak, as 
realized also in Fig. 5.2 for the lowest T = 0.1 t. 





Figure 5.2: Dynamical mobility (i{uj) for a single hole in an AFM at various T. 

From u{ui) one can extract the d.c. mobility (1q. Results for J = 0.3 t and J = are presented 
in Fig. |5.3j , with the RPA (J = 0) result for comparison. While for T > 2t we see the high-T result 
(jlq oc 1/T for both J, we get for the AFM case a nontrivial /xo ~ const within the regime J < T < t. 
It is plausible that J > reduces (j,q as follows also from the frequency-moment analysis at T> f 
(Jaklic and Prelovsek 1995c), indicating that in an AFM the hole motion is more frustrated leading 



to stronger scattering. As one can realize from Fig. 5.2, it is hard to get meaningful results for 
T < J ~ Tf s , since in this regime the central QP peak is essentially undamped. This is an indication 
that we are already dealing with finite-size effects, as discussed in Sec. 5.2, and the mean free path 
is beyond our system size l s > L. Naively one would expect from Fig. |5.3| that in the regime T < J 
(jlq would increase with decreasing T, whereby the scattering of the QP (AFM polaron) should be 
on thermally excited magnons. There are however unsolved problems with such a description. In 
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an AFM the QP dispersion (Kane at al. 1989, Martinez and Horsch 1991) is strongly renormalized 
and is effectively narrower than the magnon one, hence it seems hard to find appropriate allowed 
scattering processes. These questions become important when discussing the relation of calculated 
to the resistivity of cuprates at low doping, as discussed in Sec. 5.5. 




5.4 Optical conductivity at finite doping 

Even more challenging questions appear at finite hole doping. Whereas at low doping one could 
treat the transport within the semiconductor-like model of independent QP - AFM polarons, this 
concept fails even for a moderate Ch ^ 0.1 due to the overlap between extended spin deformations 
around holes. Alternative indications for such a phenomenon have been already discussed in Sec. 4 
in connection with Ch(nh)- 

Since in a planar system we are evaluating the sheet conductivity ct(lo) it is convenient to discuss 
at finite doping the dimensionless quantity a = crpo, where po = h/e^ = A.lkCl is the universal 2D 
sheet resistance. Such a quantity has an additional meaning since it makes a direct contact with 
the theory of localization, where ao ~ o m in is a characteristic marginal value associated with the 
2D minimum metallic conductivity (Mott and Davis 1979). The value of a m i n and its relevance is 
however controversial, and it ranges from a m i n = 0.1 (Mott and Davis 1979) to a m in ~ 0.5 found 
in experiments (Mandrus et al. 1991) as a borderline between insulating and conducting cuprates 
at low T. 

5.4.1 Intermediate doping 

Let us go straight to results at the intermediate (optimum) doping, where we can reach lowest Tf s . 
Instead of cr(ui) it is more instructive to present the current correlation function C(lo), equation 
( |5.1| ). To avoid the ambiguities with an additional smoothing, we plot the corresponding integrated 
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spectra 

l c (u) = I C{u')dw'. (5.12) 



o 

In Fig. |5,4| we present Ic{u) for = 3/16, for various T < t. Spectra reveal several remarkable 
features (Jaklic and Prelovsek 1995b): 

(i) For T < J spectra /c(^) are essentially independent of T, at least for available T > Tf s . 

(ii) Simultaneously the slope of Ic{^ < 2 t) is nearly constant, i.e. we find C{uo) ~ Co in a wide lo 
range. At the same time Co is only weakly dependent on J as tested for J/t = 0.2 — 0.6. 

(iii) Even for higher T > J the differences in the slope Co, as also in Ic(oo), appear as less essential. 
Note that for T>fwe know exactly Ic(oo) = irt 2 Ch(l — Ch) (Jaklic and Prelovsek 1995c). 

We conclude that C(uj < 2 t) ~ Co implies a simple universal form (Jaklic and Prelovsek 1995b), 

1 — e _/3w 

a(u) = C e 2 . (5.13) 

Such a c(uj) shows a nonanalytic behaviour at to — > 0, starting with a finite slope. This is already 
an indication that cr(u;) is not consistent with the usual Drude form, but rather with a marginal 
concept (Varma et al. 1989) where the only to scale is given by T. It is also remarkable that the 
form ( |5.13j ) trivially reproduces the linear law pocTas well as the non-Drude fall-off at uj > T. It 
is evident that the expression (5.13) is universal containing the only parameter Co as a prefactor. 



Experimental results and theoretical considerations are often discussed in terms of the in- 
dependent relaxation time r and the effective mass m* . These can be uniquely introduced via 
the complex d{oj) and the corresponding memory function M(u) (Gotze and Wolfle 1972), 

^ = Z^lk7y S = -(H km )/2N, (5.14) 
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Figure 5.5: (a) Inverse relaxation time 1/r, and (b) mass enhancement m* /mt vs. u, for Ch = 3/16 
and different T. The insert in a) shows the T-dependence of l/r(0). 

and 

1 M"(uj) 



t{u) 1 + M'(u)/w' 
m» _ 2c,^ i+ MV) 



where mt = \j1ta\ is the bare band mass. Using relations ( |5.14 ),( 5.15| ) one can formally rewrite 
(t(uj) in the familiar Drude form 

m*(u)[u + i/t(lu)\ 



Employing the equation ( 5.16 ) we evaluate both t(u) and m*(tu) from known Results for 

Ch = 3/16, corresponding to Fig. p.4| , are presented in Fig. |T5|. 

It follows from Fig. |5.5| a, but even more directly from the form ( |5.13| ) , that in the regime T < J 
and to < t we can describe the behaviour of 1/r with a linear law 

t- 1 = 2tt\(lo + (,T), A ~ 0.09, ^~2.7. (5.17) 

This dependence falls within the general framework of the MFL scenario. It is however not the 
form ( |2.3D proposed originally (Varma et al. 1989, Littlewood and Varma 1990), but rather the 
one (|2.4D deduced from experiments in cuprates (El Azrak et al. 1995, Baraduc et al. 1995). It 
should be however stressed that the asymptotic form ( |5.17 ) does not allow for any free parameter, 



i.e. constants A and £ are universal and independent of any model parameters, whereas within the 
MFL proposal A is an adjustable parameter while £ = ir. 

It should be also mentioned that the universal dynamics, as described by the equation (|5.13j ), 
does not seem to be restricted only to the particular case of doped AFM, but has a wider ap- 
plicability, e.g. it has been recently established also in ladder systems (Tsunetsugu and Imada 
1997). 
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5.4.2 Underdoped and overdoped regime 



Let us turn to the discussion of results at other dopings c^. Results for = 4/16 (Jaklic and 
Prelovsek 1995c) are in all aspects very similar to the Ch = 3/16 case. We show in Figs. |5.6| analogous 
Ic(u) for underdoped = 2/18, as well as for overdoped = 7/16. One feature common to all 
Ch considered (including = 1) is the non-Drude behaviour for u > t. This confirms the belief 
that the incoherent motion of holes, dominating the high-w response as described e.g. within the 
RPA (Rice and Zhang 1989), can remain a valid concept even at large doping Ch < 0.5. 




Figure 5.6: Ic(w) for (a) Ch = 2/18, and (b) Ch = 7/16, for different T. 
Otherwise it is harder to make conclusions on the most interesting low-(u;,T) regime, both for 



the underdoped as well as for the overdoped case. For = 2/18 we note that in contrast to Fig. 5.4 
the low-o; slope is not constant at T < J but is still gradually decreasing with T. This results in 
a modified variation of the d.c. conductivity o"o oc l/T 1 ^ with r\ > 0, hence also in a sublinear 
po oc T 1 "^ . On the other hand we see already for T = 0.1 t a drop of Ic{oo) (closely related to 
the sum rule) which can be interpreted with a significant persistent (nonscattered) current in a 
system with p.b.c. The latter emerges as the coherent part with D c > according to the equation 
( |5.3| ), not taken into account in the presented Ic(oj). While we cannot say much on the current 
relaxation of such a coherent part, it is plausible that an increase in the coherence should decrease 
the resistivity po, as found in underdoped cuprates below the crossover T*. Otherwise results in 



Fig. 5.6 a are close to curves for a single hole in an AFM, Fig. |5^, with the difference that in the 
latter case T* appears even higher. 

In the overdoped case in Fig. |5.6| b Ic{oo) starts to deviate from a constant already for T < J. 
At T < J a sharp increase of an unscattered current appears, so we are unable to speculate on a 
possible onset of more LFL-like p oc T 2 . 

5.4.3 Effects of the next-nearest-neighbour hopping 

It is a relevant question whether the anomalous but universal behaviour of cr(u>) found at the in- 
termediate doping is spoiled by possible additional terms, e.g. by the introduction of the n.n.n. 
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hopping t' term ( |2,7| ), invoked often to obtain a realistic description of cuprates. When considering 
the effect of t' it is important to realize that at finite doping t' > tends to stabilize the AFM or- 
dering while t' < destabilizes it and tends toward a Nagaoka-type FM state (Bonca and Prelovsek 
1989). 

Results for Ic{u) at t' / are presented in Fig. 5/7. We choose rather modest \t' /t\ = 0.2, still 
quite a pronounced effect on Ic(oo) is evident. We should note that in this case also generalized 
sum rules, following equations (|5.5[) and (|5.8[) , apply. On the other hand, the low-o; part is not 
changed essentially. For t'/t = 0.2 we note again an universal behaviour according to the form 
( 5,13p . Deviations at lowest T are somewhat larger for t'/t = —0.2. A possible interpretation of 
these results is that t' ^ moves the system effectively away from the starting doping regime, i.e. 
drives an optimum system either towards the underdoped or to the overdoped regime. 




Figure 5.7: Ic{u) for Ch = 3/18 and the n.n.n. hopping: a) t'/t = 0.2, and b) t'/t = —0.2, for 
different T. 



5.5 Resistivity 

Finally, we display results for the d.c. resistivities p(T), as extracted from a(to — » 0). The extrap- 
olation u —* is straightforward at higher T. It becomes somewhat more delicate on approaching 
T ~ Tf s due to the appearance of more pronounced structures and D c ^ 0. Actually, we can 
evaluate gq either from cr(u) which involves some smoothing, or from the slope of Ic{w ~ 0) which 
seems to be more reliable. 



In Fig. 5.8 we present results for p(T) at various Ch for J/t = 0.3, and in comparison also 
for J = 0. We first notice that at T > t calculated p{T) oc T. The slope at T > t is nearly 
independent of J, the main effect of finite J > being the upward shift of p(T) curves. The shift 
also decreases as doping is increased, since it is plausible that the effect of J nearly vanishes in 
the overdoped regime Ch > 0.3. The slope at T > t can be approximately given (for Ch < 0.25) as 
dp/dT ~ Qpoks/cht with £ ~ 0.4. This value confirms the relevance ot the RPA (Brinkman and 
Rice 1970), which yields for a square lattice Q ~ 0.72. 
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Resistivity p vs. T at different Ch = Nh/N for J = 0.3 t (full lines) and J = (dashed 



One could expect essential changes in the regime T < J. It is evident from Fig. |5.8| that at the 
intermediate doping 0.15 < c\ t < 0.25 the curve p(T) changes the slope at T ~ J. Nevertheless 
p(T < J) is still linear in T, as emerges also from the observed universal form ( |5.13| ). As discussed 
in Sec. 5.4.2, the underdoped systems as Ch = 2/16 show already deviations from the universality 
( 5. 13] ). Two effects should be mentioned. In the regime where our analysis is fully applicable and 
l s < L we find sublinear p(T) oc T 1_r? , which could be interpreted also as p(T) ~ A + BT with a 
positive intercept A > 0. On the other hand, the appearance of D c ^ in our results seems to be 
closely related to another scale T ~ T* where a kink of p{T) appears. It should be also noted that 
the introduction of the n.n.n. hopping t 1 ^ does not change qualitative conclusions. 



5.6 Relation to experiments 

In 2D a is naturally expressed in terms of the universal constant po = fr/e^. The corresponding 
3D conductivity of a stack of 2D conducting sheets with an average distance c is given instead 
by 03D = cr/c. For comparison with experiments we reduce the 3D measured values to the 2D 
conductivities for three different cuprates at intermediate doping, i.e. LSCO with x = 0.2 (Uchida 
et al. 1991), Bi 2 Sr 2 CaCu 2 8 (BISCCO) (Romero et al. 1992), and YBCO (Cooper el al. 1993). 
Experimental results are taken at relatively low temperatures T < 200 K, i.e. at T < Tf s below the 
sensitivity of our small-system calculations. It is an open experimental question whether spectra 
reduced in this way are really universal, i.e. whether they depend only on the hole concentration 
Ch in conducting Cu0 2 layers or other details of material structure are still important. Moreover, 
effective for presented materials are known only approximately, except for LSCO where ~ x, 
i.e. it is estimated that c h ~ 0.23 for BISCCO and YBCO (Batlogg et al. 1994). 

Taking into account this uncertainty calculated spectra cr(u) for = 3/16 are in a quantitative 
agreement with measurements, as seen in Fig. 5.S, where the energies are expressed in eV using 
t = 0.4 eV. We note also that at high uj > 1 eV the fall-off of calculated (j{uj) is faster than that 
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of measured ones. This could be explained by the emerging contribution of transitions to excited 
(charge transfer or upper Hubbard band) states not taken into account within the t — J model, but 
clearly identified experimentally (Uchida et al. 1991). 



n h =3/16 
J/t=0.3 
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Figure 5.9: Sheet conductivities (t{uj) for various T/t, in comparison with measurements in different 
cuprates. Experimental results refer to T < 200 K. 

It should be noted that also calculated r _1 (w) and corresponding parameters A and £, as defined 



by the equation ( 5.17 ) are close to the experimental ones. E.g., an analysis in a wide frequency 



range u < 1000 cm" 1 presented for YBCO and BISCCO data by El Azrak et al. (1994) yields 
A ~ 0.11, while a broader class of optimally doped materials is consistent with 0.1 < A < 0.15 
(Baraduc et al. 1995). A recent investigation (Startseva et al. 1997) in a narrower to range obtains 
for an overdoped LSCO (x = 0.22) similar values 0.06 < A < 0.08. 

Also m* /nit in Fig. [5l]b qualitatively agrees with experimental findings (Romero et al. 1992, 
El Azrak et al. 1994, Tanner and Timusk 1992, Cooper et al. 1993). For a quantitative comparison 
we should note that by lowering T to experimentally investigated values we expect the increase of 
m* at low uj. On the other hand, assuming realistic values ao = 0.38 nm and t = 0.4 eV we should 
also take into account that nit ~ 0.6 niQ. 



In Fig. |5.10| we compare calculated resistivities to the measured ones. It should be pointed out 
that there is a restricted T window where a comparison could be made since Tf s ~ 450-ff , whereby 
at T ~ Tf s also finite-size effects start to influence our analysis. Nevertheless, for the intermediate 
doping Ch ~ 0.2 our p(T) results match quantitatively well experimental ones for cuprates with 
comparable hole concentrations, i.e. for BISCCO, YBCO and LSCO with x = 0.15 (Takagi et al. 
1992, Batlogg et al. 1994, lye 1992). 

Turning to underdoped materials, there are certain features which are reproduced in our results. 
First, our for a single hole, see Fig. |5.2| , as well as several T = numerical studies (Sega 

and Prelovsek 1990, Dagotto 1994) indicate an existence of a mid-infrared peak in ct(lo), being a 
signature of magnon excitations in a spin background with a longer range AFM order. This feature 
has been seen in experiments in LSCO (Uchida et al. 1991 ), although its appearance in other 
materials is controversial. We also reproduce the observation that at moderate doping Ch > 0.1 the 



43 



3 




1 ■ 1 ■ 1 ■ 1 

200 400 600 800 

T(K) 

Figure 5.10: Sheet resistivities p(T) for various dopings (full lines) in comparison with measure- 
ments in LSCO with x = 0.15 (dotted), BISCCO (dashed), and YBCO (dash-dotted). 



resistivity essentially scales as p oc 1/c^. An unproportional increase of p at low doping < 0.1, 
as deduced when comparing values in Fig. 5.1C with single-hole mobilities in Fig. 5.3, is as well 
consistent with experiments, e.g. for LSCO with x < 0.1 (Takagi et al. 1992). 



6 Magnetic Properties 

The static spin response and the spin dynamics in the undoped and doped AFM state in cuprates 
have been experimentally studied by the neutron scattering (see e.g. Shirane 1991), by NMR 
techniques (see e.g. Slichter 1994), as well by other methods. 

Magnetic excitations in cuprates as measured by the neutron scattering reveal in insulating 
materials La2Cu04 (Keimer et al. 1992, Hayden et al. 1996) and YBa2Cu30g (Rossat-Mignot 
et al. 1991) a remarkable agreement with the magnons within the planar Heisenberg AFM. On 
the other hand, a consistent description of magnetic properties of doped cuprates is still lacking. 
Nevertheless it is quite clear that the normal-state spin dynamics differs qualitatively from the one 
expected for LFL. In doped LSCO NMR and NQR spin-lattice relaxation time T\ on Cu nuclei 
is generically nearly T and doping independent in the normal state (T > T c ) (Imai et al. 1993), 
in contrast to the Korringa law \/T\ oc T in normal metals. Also, in the same regime the low-u; 
dynamical susceptibility in doped systems appears to be consistent with x"( w ) oc w/T (Shirane 
1991, Keimer et al. 1992, Sternlieb et al. 1993). Experimental investigations in recent years have 
been focused on underdoped materials, which show in the normal phase the emerging spin gap 
both in the NMR (Slichter 1994) and in the neutron scattering (Sternlieb et al. 1993), but also the 
appearance of charge stripes related to AFM domains (Tranquada et al. 1995). Since the latter 
phenomena appear only at lower T, hardly accessible in our studies, we shall comment them only 
briefly later on. 

Well understood theoretically is so far the isotropic quantum AFM, with the long-range order at 
T = (Manousakis 1991). For doped systems several phenomenological explanations (see Sec. 2.2) 
have been presented for magnetic properties. NMR and NQR data on the spin dynamics have been 
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interpreted within the NAFL scenario (Millis et al. 1990, Millis and Monien 1992), where the T 
dependence is attributed to the variation of the AFM correlation length £(T). At low hole doping 
the mapping on the quantum critical scaling regime of the nonlinear sigma model, with the main 
ingredient £ oc 1/T, has been advocated by Sokol and Pines (1993). An alternative scenario for the 
1ow-(cj,T) behaviour is the MFL scenario (Varma et al. 1989, Littlewood and Varma 1991), where 
£(T) is short and not crucially T-dependent. 

Most reliable model results for magnetic properties of doped AFM have been obtained for static 
quantities by QMC studies of the Hubbard model and by the g.s. ED of the t-J model (see Dagotto 
1994), as well as by means of the HTE analysis (Singh and Glenister 1992a). Much less conclusive 
are results for the spin dynamics. While g.s. ED studies of the spin structure factor S(q,tv) 
(Tohyama et al. 1995, Eder et al. 1995) reveal quite an essential difference between S(q,uj) and 
corresponding charge spectra N(q, u), they cannot give reliable conclusions for the low-w behaviour. 
There have been only few attempts to address numerically dynamical properties at T > (Tohyama 
et al. 1993), prior to the application of the FTLM (Jaklic and Prelovsek 1995a, b). It should be 
reminded that most challenging questions (related to the normal state) refer to the dynamics at 
low ijj and to the d.c. spin response in the strong correlation regime J < t at low T < J. 



6.1 Spin response 

We consider the response of the electronic system to the time-dependent, spatially modulated 
magnetic field, which couples to the spin degrees of freedom via a Zeeman term 

i 

Within the linear response theory the magnetization response is given by 

<y(M|)(w) = g 2 n 2 BX (q,u)B g iu;), (6.2) 

where 

POO 

X (q,u) = i / dte^([SI(t),S^(0)]), 
Jo 

SI = (1/VAO ^V f %f. (6.3) 

i 

x"(q,w) can be expressed with the dynamical spin structure factor S(q,uj) 

X "(q,cj) = 7r(l-e-^)S(M, 
1 f°° 

S(q,u) = -Rej^ dte^(SI(t)S z _^0))- (6-4) 

S(q,u) is the quantity directly evaluated within the FTLM, following the expression ( 3.42j ), Con- 
sequently the static structure factor and the static susceptibility can be evaluated, 



S(q) = / S(q,u;)du; = (SiS^), 

J —oo 

X(q) = x!im = -T > r^^-du. (6.5) 



71" ./-oo <*> 
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6.2 Uniform spin susceptibility and Wilson ratio 

The uniform case q = is particular, since 5"*°* is a conserved quantity. Then one can use for the 
uniform static spin-susceptibility xo the expression 

/(Qtot\2\ 

Xo = X(Q = 0) = PS(q = 0) = (6.6) 

The calculation reduces to the one involving only a thermodynamic average of a conserved operator, 
and an analysis analogous to the one used in Sec. 4 applies. This simplification allows for the 
consideration of larger systems, i.e. we calculate Xo{T) for a system of N = 20 sites in the range 
< c/j < 0.3, while for the undoped AFM (Heisenberg model) we reach N = 26 (Jaklic and 
Prelovsek 1996). 



Results at several are presented in Fig. 3.1. For = our results agree with the HTE (Singh 
and Glenister 1992a) down to T ~ 0.3 J. In an AFM xo exhibits a maximum at T = T* ~ J, which 
reflects the onset of the short range AFM order for T < T*. Namely, at T — > the longitudinal 
spin fluctuations in an ordered AFM gradually freeze out, while transverse ones remain constant, 
thus leading to a reduced Xo(T — > 0). 

For a finite doping xo is only weakly diminished for T > J, the reduction being due to a reduced 
number of spins in the system. On the other hand, the qualitative change appears for T < J. The 
maximum T* gradually shifts to lower T with doping and finally disappears at Ch > 0.15. In the 
overdoped regime Ch > 0.15 we observe a monotonous Pauli-like xo f° r T < 0.2 t, which could 
signify an onset of a low-T behaviour consistent with the LFL picture. Still up to c/j > 0.6 we do 
not find the usual LFL behaviour Xo(^) = const, but rather x(T) still increases nearly linearly on 
lowering T — > 0. 

Obtained results are quite consistent with experiments on cuprates. First, xo a t T — > increases 
with doping as found in LSCO (Johnston 1989, Torrance et al. 1989) and the existence of a 
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Figure 6.2: Wilson ratio Rw vs. T at several Ch- 



maximum of Xo(^) at T*(ch) in underdoped cuprates has been interpreted in terms of a pseudogap 
scale (Batlogg et al. 1994). It is evident from our results that such a pseudogap is inherent within 
the t-J model and is intimately related to the short range AFM order. It is remarkable that in 
the model this feature disappears, i.e. T* — > 0, just at a similar (optimum) doping ~ 0.15 as 
in cuprates (Batlogg et al. 1994). Note also that the unusual linear increase of Xo(T — > 0) in the 
overdoped regime, see Fig. [O] for Ch = 0.2, seems to be well consistent with experimental results 
in overdoped LSCO (Loram et al. 1996). 

Here we can also comment on the ratio W = (g 2 ^ 2 B Xo)/l where 7 = Cy /T. It is used as a test 
for the concept of nearly free QP, where one expects Wq = ^(irks / ' Hb) 2 • The meaningful measure 
is thus the Wilson ratio Rw = W/Wq being often studied in correlated systems, in particular in 
connection with heavy-fermion metals. In our case it is convenient to perform the test on s directly 
by defining 7 = s/T (see also Loram et al. 1996). In our notation the dimensionless Wilson ratio 
can be expressed as 

4vr 2 (k B T/t)( X ot) 



Rw 



s/k 



B 



(6.7) 



We can now easily evaluate Rw by comparing results in Fig. 4.3a with Fig. 6.1 and results are 
shown in Fig. |6.2| . It is quite striking that in the intermediate regime 0.1 < < 0.2 at low T <^ J 
we find values, very close to the free fermion one, i.e. Rw = 1- In the overdoped regime calculated 
values are somewhat larger, but still Rw < 2. The same seems to be the case for the underdoped 
case Ch = 0.05, where results at lowest T should be taken with care. With increasing T also Rw is 
increasing and seems to reach quite a wide plateau for T > J with Rw ~ 2. 

This finding is quite puzzling since clearly doped AFM are far from a simple LFL, and apparently 
even further from free fermions. Note however that experimental results also yield Rw ~ 1 ; both 
for LSCO and YBCO (Loram et al. 1996). Moreover, the same experiments indicate that s cx Txo 
in a wide range of T and doping, so Rw as defined in the equation Q6.7| ) is nearly T independent. 
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6.3 Spin structure factor and dynamical susceptibility 

Let us first discuss S(Q,u) spectra at the AFM wavevector Q = (vr,7r), where the spin response 
is most pronounced. In Fig. |6.3| we present results at fixed T = 0.2 t < J, but various dopings 
c/i (Jaklic and Prelovsek 1995a). It should be noted that S(Q,oj) is not symmetric around oj = 0, 



hence a maximum appears in general at oj > 0. The most interesting feature in Fig. |6.3| is the 
qualitative change of spectra on doping. At low doping < 0.12 we see that S(Q, oj) is dominated 
by a single central peak, i.e. we are observing mainly the AFM spin fluctuations of a Heisenberg 
model. The main effect of holes is to reduce the AFM correlation length, and consequently the 
intensity of the central peak. 

In the intermediate regime 0.12 < Ch < 0.3 a high-frequency component with oj > t emerges, 
coexisting with the remaining low-o; fluctuations. It is quite plausible to attribute the high-tj 
dynamics to the free-fermion-like component of the correlated system, in particular since it appears 
to be quite independent of J (provided that J < t). Although this observation is consistent with 
previous studies (Tohyama et al. 1993, Putikka et al. 1994), the coexistence of spin-fluctuation and 
free-fermion timescales at the intermediate doping has been clearly established only by using the 
FTLM (Jaklic and Prelovsek 1995a). Namely, it is harder for other methods, e.g. within the HTE 
method, to resolve coexisting different timescales. The dual character is a crucial property, since 
the free-fermion part determines to large extent static spin correlations S(q) through the relation 
(|6.5|), as well as the electron (charge) density correlations N(q), discussed in Sec. 8.1. The latter 
have been in fact interpreted in terms of a quasi FS (Putikka et al. 1994). On the other hand the 
low-u; spin dynamics dominates dynamical and static spin susceptibilities, i.e. x" {q,oj) /oj and x(,Q)i 
hence the low-w neutron scattering and NMR processes. It is also evident that in the overdoped 
regime Ch > 0.3 the low-frequency component is disappearing and the free-fermion-like fluctuations 
tend to exhaust the spectra. S(Q,oj) is nevertheless quite distinct from the one for free fermions, 
even at c/j ~ 0.5. 

Before we present results for general q ^ Q, let us discuss the AFM correlation length £(T). 
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One can evaluate £(T) from the static real space correlations S(r), corresponding to S(q), i.e. 

^ = 7^ E \^ ex PW • ( 6 - 8 ) 
45(Q) Y 

In the most interesting regime Ch = 0.1 — 0.3 we find that £ is short, typically £ ~ 1, governed by 
correlations at r% = 1. It increases less than 30% between T = J and T = J/3. This finding (Jaklic 
and Prelovsek 1995a) is well in agreement with the HTE studies (Singh and Glenister 1992a), with 
the QMC results for the Hubbard model (Furukawa and Imada 1992) as well as with values £(T = 0) 
obtained via the ED within the t — J model (Dagotto 1994). Similar values for £ can be extracted 
also considering our results for static x{<t)i although results are less conclusive at low T ~ Tj s . 

In recent years experiments in cuprates have shown the possibility of longer range (or even 
long range) spin correlations with q ^ Q (Tranquada et al. 1995, Hayden et al. 1996). This has 
been related either to an incommensurate spin order or to the appearance of charge stripes. Such 
structures are well possible within the t-J model (Prelovsek and Zotos 1990, White and Scalapino 
1997b) although it is controversial whether they are stable in the most interesting regime J <C t. 
Nevertheless we do not expect to resolve them in our calculations for J/t = 0.3, since they could 
possibly appear only for T < Tf s and could be also missed due to particular clusters which do not 
prefer expected orderings. 

Let us turn to dynamical susceptibilities. Fig. |6^J displays \"(q,uj)/oj for fixed Ch = 3/16, but 
various T and nonequivalent q. Note that on a 4 x 4 lattice q = (0,7r) and q = (ir/2, n/2) are 
equivalent. In contrast to Fig. |6.3| , high-u; features are now suppressed. Also, the free-fermion-like 
component is well separated from the 1ow-cj part only for q ~ Q, where one expects a gap in the 
response of free fermions with a well defined FS. For other q the free-fermion contribution persists 
at larger to > J in the form of a long tail, while in the low-o; regime it merges with the spin 
contribution. 

The most striking feature of Fig. |6.4| is the strong T dependence of the low-o; spectra, whereas 
the AFM correlation length £ has been found to be only weakly T dependent. This conclusion on 
x"(q,u>)/u! seems to hold within the correlation regime \q — Q\ < £ , where also x(o) ~ x{Q)- The 
relevant volume in the q space clearly increases on doping and exhausts for = 3/16 already the 
substantial part of the Brillouin zone. The variation at low u appears to follow, 

X"(q,u)/u oc 1/T, uj<T<J, (6.9) 

or equivalently, from the relation Q6.4j ) S(q,uj) is nearly T and uj independent in the same regime. 
On the other hand, at the same doping scaling does not hold for q outside the correlation volume, 
e.g. for q = (0,7r/2) in Fig. 6.4, where x"(<f>^)/^ is approximately T independent as expected 
within the LFL. 

Spectra discussed above have as a direct consequence the T variation of the static x(o) f° r 
T < J. We observe a pronounced T dependence, e.g. x(o) ^ 1/T in a wide regime J/3 < T < t 
for all q within the correlation regime. It should be however noted that we are quite restricted in 
the range of T/J, so that more quantitative conclusions on a possible power-law (or logarithmic) 
variation with T are not feasible. 
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Figure 6.4: x" {Q-^^/u) a t c h = 3/16 for nonequivalent q and different T: T/t = 0.1 (full line), 0.2 
(dashed line), 0.3 (dash-dotted line), and 0.5 (dotted line). 



6.4 Local spin dynamics 

In order to describe properly the spin correlation function S(q,u>) it is very helpful to consider the 
local spin correlation function Sl(uj) and its symmetric part S(uj) (Jaklic and Prelovsek 1995b), 

q 

S(u) = S L (u;) + S L (-u;) = (l + e-^)S L (u;). (6.10) 

It should be noted that Sl(u) and the related susceptibility xl( w ) have been directly measured in 
cuprates by the neutron scattering (Shirane 1991), while the NMR relaxation on 63 Cu effectively 
yields the information on Sl(w — ► 0), as discussed in Sec. 6.5. An important restriction on S(u) is 
the sum rule following from the equation (|6.5| ), 

OO 1 

S(u;) ( hj = ((S!) 2 ) = -(l-c h ). (6.11) 

When we perform the calculation of S(lo) we omit the singular q = term in the summation ( |6.10| ), 
so the sum rule (|6~lll) serves as an useful test. In Fig. pJ3 we display S(u) for different dopings 
Ch = 1/20 and 3/16, and several T in the range 0.1 < T/t < 0.7 (Jaklic and Prelovsek 1995b). It is 
immediately evident that S(u>) at the optimum doping Ch = 3/16 is essentially T independent in a 
wide T range, although one crosses the exchange-energy scale T ~ J. 
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Figure 6.5: Local spin correlation function S(u>) for = 1/20 and = 3/16, and different T. 



For the underdoped case = 1/20 as well as for the undoped AFM the behaviour is analogous 
for higher T > To ~ 0.7 J. Deviations appear at T < To, leading at T —* to a decrease and to a 
flattening of S(u < 2 J) ~ const, whereby a weak enhancement of S(uj > 2 J) is then required by the 
sum rule. At least this regime of essentially pure Heisenberg model is well understood theoretically. 
The behaviour at T > To is consistent with the quantum critical regime within the AFM, whereas 
for T < To the renormalized classical regime applies (Chakravarty et al. 1989). In the latter regime 
we are dealing with longer range AFM correlations £ 3> 1, hence S(oj) is essentially that of an 
ordered AFM where the simple magnon picture leads in 2D to S(uj < 2 J) ~ const. 

In Fig. 3^ we follow the doping dependence of S{u>) at fixed T = 0.2 t < J. For convenience 
we plot again the integrated spectra, in analogy to the equation ( 5.12j ), 



S(oj')duj' . 



(6.12) 



For chosen T results are again most reliable at the intermediate doping. The most striking message 
is that the initial slope of Is(w) and consequently Sl{uj — > 0) is nearly doping independent for 
< Ch < 0.25, as well as T independent at the intermediate doping. Only in the overdoped systems 
with Ch > 0.25 the low-w behaviour changes qualitatively, where the low-w contribution is strongly 
suppressed as expected in a (more) normal FL. In a pure Heisenberg model the spin dynamics is 
nearly exhausted in the range u < 3J with excitations being magnons. On the other hand, at 
the intermediate doping S(u>) decreases smoothly up to u ~ At, this being a consequence of the 
free-fermion-like component. Still up to Ch ~ 0.3 the dominant scale of spin fluctuations remains 
related to J. 

Using equations ([T3) and ( |6.10| ) we finally note that the T-independent S(u>) leads to the 
universal form 

x '>) = itanh(|^)S(u/), (6.13) 
where the T dependence is only in the thermodynamic prefactor. It follows from the equation 
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Figure 6.6: Integrated spin correlation spectra Is(w) at fixed T = 0.2 £ and different = Nh/N. 



( 6.13| ) that for T < J the only relevant scale for Xz,( w ) is given by T. The same should hold for 
the response at fixed q. So one can generalize the expression ( |6.13D to 

27rln" 1 (£g m ) 



XL 



k--Qi 2 +r 2 



(6.14) 



where the cutoff g m ~ 7T. The scaling is expected to be valid at least for \q — Q\ < Since 
£(T) remains weakly T dependent on decreasing T, this introduces an additional T variation in 
the expression ( |6.14 ). Far from the regime q ~ Q, in particular for q ~ 0, the response is more 
free-fermion like, i.e. x"{q^) 1S T independent, as seen in Fig. 6.4 for q = (tt/2, 0). 



6.5 Nuclear magnetic relaxation 

For experimental studies of static and dynamical spin correlations the NMR and the NQR relaxation 
are among the most valuable tools. The hyperfine coupling of electronic spins to 63 Cu and 17 
nuclear spins / within the CuO plane has been established by a number of authors (Mila and Rice 
1989, Shastry 1989, Millis et al. 1990, Millis and Monien 1992, Slichter 1994), 



(6.15) 



Form factors are determined by the relative position of electronic spins (assumed to be centered on 
Cu sites) and nuclear spins, 



63 A(<z) 
17 A a (q) 



A + 2B(cos q x a + cos q y a), 
2C cos(q a a/2), 



(6.16) 



where a = x,y refer to Cu-O-Cu axis orientation, and A, B, C are the hyperfine couplings. A is 
the direct hyperfine tensor, coupling the electronic and the nuclear spin on the same Cu site, with 
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components A±_ and Ah, corresponding to spins oriented perpendicular and parallel to the Cu02 
plane, respectively. 

The nuclear spin-lattice relaxation is due to the coupling ( 6.15| ), directly related to the low-w 



electronic spin fluctuations. Let us consider only the magnetic field directed perpendicular to the 
CuC>2 plane. Then the NMR spin-lattice relaxation rate is given by (Slichter 1994), 

^ = 2^El^®| 25 ®^o), (6-17) 
q 

where loq is the precession frequency of the nuclear spin in the magnetic field, and £ = 1. Since 
ujq -C T, one can express result as well with the dynamical susceptibility S(q,u>) Tx"(q, ujq)/ujq. 
The NQR relaxation rate 1/Tl is given also by the expression fl6.17| ), but with £ = 4 (Millis and 
Monien 1992). Form factors ( |6.16| ) differ essentially for O and Cu sites. Due to AFM fluctuations 
spin fluctuations are at maximum around q = Q, as in the equation ( |6.14 ). The weight of this region 



enhances 1/Xi for 63 Cu, while suppressing it for 17 0. This is the most important point explaining 
the observation of very different relaxation rates on different nuclei. 

In the evaluation of 1/Ti within a finite system we again omit the q = term since x"(q> w ) 1 00 i s 
ill defined for u — > due to conserved total S z . A proper treatment would require a more detailed 
spin-diffusion contribution at q ~ which however seems to be less important, at least for the 
undoped system (Sokol et al. 1993). To allow a direct comparison with experiments we choose 
hyperfine-coupling parameters as proposed in the literature (Millis and Monien 1992). 

Results for 63 Cu NQR relaxation rate 1/Ti are presented in Fig. |6.7| (Jaklic and Prelovsek 
1995a). For the undoped case our results agree with Sokol et al. (1993). It is remarkable that 
1/Ti appears to be nearly T independent for a broad range of hole concentrations 0.06 < Ch < 0.3. 
Only for overdoped systems with > 0.5 the behaviour at T < t approaches that of a normal LFL 
with 1/Ti of T. Since the 63 Cu form factor is maximum (but slowly varying with q) near AFM Q, 
we get approximately 1/Ti °c Sl(wo ~ 0). Previously established universality of Sl(u>) naturally 
explains nearly T independent 1/Ti, which is moreover also weakly dependent on Ch for Ch < 0.15. 

Results as shown in Fig. [T?] are in agreement, even quantitatively without any fitting parame- 
ters, with remarkable NQR experiments on doped LSCO (Imai et al. 1993), which reveal nearly T- 
and x-independent 1/Ti at T > 300 K and x < 0.15. In contrast to experiments our results show 
some variation of 1/Ti with Ch in this doping range. For the optimum doping lower rates 1/Ti are 
anyhow expected, consistent with the data for YBCO (Takigawa et al. 1991, Kitaoka et al. 1991), 
where for T > T c the rate 1/Ti is again only weakly T dependent. 

Essentially different T dependence of 1/Ti on Cu and O nuclei has been used as an evidence for 
the importance of strong AFM correlations and for the non-LFL behaviour in doped cuprates. To 
evaluate the NMR relaxation rate 1/Ti for 17 we can use again the expression (fTl7D with proper 



form factors (6.16), projecting out the AFM fluctuations at q ~ Q. The omitted q ~ contribution 



introduces in this case a larger uncertainty. Nevertheless, for Ch = 1/16 and Ch = 2/16 we recover 



results very well described with the Korringa behaviour, i.e. 1/Ti ~ wT, as seen in Fig. 3. J for 
Tf s < T < J. In particular for Ch = 2/16 we get w ~ 0.3(sK) -1 , very close to the actual value 
w ~ 0.35(sK)~ 1 reported for YBCO (Takigawa et al. 1991, Kitaoka et al. 1991). For c h > 3/16 
deviations from the Korringa law become more pronounced due to very short £ ~ 1. 

The Gaussian component of the spin-spin relaxation rate I/T2 for 63 Cu nuclei can be on the 
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Figure 6.7: 63 Cu NQR spin-lattice relaxation rate 1/T\ vs. T for different dopings c^. 




other hand related to static spin susceptibilities (Slichter 1994), 



T, 



, 1/2 

' TV 



(6.18) 



The ratio R = T\T '/T2 for 63 Cu nuclei relaxation, which is approximately T independent in cuprates, 
has been interpreted as the evidence for the quantum critical behaviour of the effective spin system 
(Sokol and Pines 1993). We find quite analogous T variation of R (Jaklic and Prelovsek 1995a) 
using calculated x(<f, u) and equations fl6.17|) , ( |6.18 ), together with the standard parameters (Millis 



and Monien 1992). The origin of R(T) ~ const is however considerably different from the quantum 
critical scenario, since the anomalous T2(T) dependence can be related to the T dependence of 
static x(q) which does not seem to be connected in an evident way with the variation of £(T). 
Results indicate stronger doping dependence, even at low doping. Quantitatively, obtained values 
are in reasonable agreement with the experimental ones, e.g. R ~ 1700 K at T = 300 K for YBCO 
(c h ~ 0.23), while R ~ 2400 K for YBa2Cu 3 6 . 63 (Sokol and Pines 1993). 

6.6 Neutron scattering 

The standard neutron scattering using the thermalized neutrons from the reactor is restricted to 
investigations in the energy regime oj < 40 meV. For cuprates this means only the low-energy part 
of the spin dynamics with u < J. There have been several measurements, probing the local Sl(co) 
by measuring S(q,u) integrated over q ~ Q (see Shirane 1991). To compare with our results one 
can simplify the relevant expression fl6.13[ ) by S(u>) ~ Sq. Such a form has indeed been used to 
describe experiments (Keimer et al. 1992, Sternlieb et al. 1993). In this connection we should 
take into account that we are not able to establish in our model calculations (the main reason 
being too high Tf s ) the existence of the pseudogap w 9 ~ 0.1 J, observed in cuprates at low T >T C 
(Rossat-Mignod et al. 1991, Sternlieb et al. 1993) and invoked in the detailed form of Sl(u>) used 
in describing neutron scattering experiments. 

The introduction of spallation neutron sources has expanded crucially the accessible energy 
range of the spin dynamics to u < leV. This gives a possibility to measure the whole relevant 
range of the spin dynamics in cuprates, results being available for undoped La2Cu04 and SC 
LSCO with x = 0.14 (Hayden et al. 1996). Data for the undoped material are well described with 
the standard magnon excitations in an ordered AFM, so more challenging are results for the doped 
material. Since the authors present the calibrated local x'li 1 - ) up to w ~ 0.25 eV, we can compare 
data quantitatively with our results e.g. for Ch = 3/16, as shown in Fig. [3j| For the latter case 



measurements were performed at T = 17 K<C J, so we present in Fig. 6.9 the T — > limit of our 



result following equation ( 6,13j ) x'l^ > T) ~ S{uj)/it. We observe that the agreement is quite 



satisfactory (note that we use both units in eV assuming again t = 0.4 eV) for u > 40 meV. There 
is however an additional peak at u ~ 20 meV, which does not appear in our analysis. This feature 
can be interpreted as a low-T phenomenon related to the onset of longer range incommensurate 
spin order found in the same material (Hayden et al. 1996), clearly beyond the reach of our method. 

6.7 Orbital diamagnetism 

In comparison with the spin response discussed in previous subsections, the diamagnetic contribu- 
tion Xd to the d.c. susceptibility has been much less investigated, both experimentally (Walstedt 
et al. 1992, Miljak et al. 1993) and theoretically (Rojo et al. 1993). A diamagnetic response 
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Figure 6.9: Local susceptibility Xzi^) as calculated for Ch = 3/16 and T 
neutron scattering result for LSCO with x = 0.14 (Hayden et al. 1996). 



0, compared with the 



emerges from the orbital motion of mobile carriers. For independent electrons it corresponds to 
the Landau diamagnetism, which is essentially T independent. For strongly correlated electrons 
such a behaviour is far from evident. Moreover, Rojo et al. (1993) have shown a relation of the 
off-diagonal Hall conductivity in an external homogeneous magnetic field B (perpendicular to the 
plane) to the orbital diamagnetic susceptibility 



a 



B 



dxd dc h 
dc h dfi ' 



(6.19) 



Since the Hall effect is known to be anomalous in cuprates, in particular its T dependence (Ong 
1990), one could speculate on similar anomalies in Xd(T), however so far both experimental and 
theoretical answers are lacking. 

Homogeneous perpendicular field B can be introduced into a tight binding model with a Peierls 
construction analogous to the equation (|5.4|), tij — > ty-exp(i0y) , where in the Landau gauge we can 
write 

% = -eoA(Ri) ■ Rij, A = B(0,x,0). (6.20) 
Our aim is to evaluate the d.c. orbital susceptibility Xd via the free energy density T 



Xd 



d 2 T 



dB 2 



(6.21) 



B=0 



It is nontrivial to incorporate in a small (tilted) lattice the phases ( |6.20 ) arising from B ^ 0, being 
at the same time compatible with the p.b.c. (Veberic et al. 1998). This can be achieved only for 
quantized fields B = tjiBq where the smallest field Bq = (j>o/N corresponds to a unit flux </>o per 
2D system. Discrete B make the calculation of Xd via equation ( |6.21 ) less reliable. It is also a 
general observation that properties involving finite B (also the Hall effect) are much more sensitive 



to finite-size effects, while at the same time translational symmetry is lost due to phases (6.20) 
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Figure 6.10: Diamagnetic susceptibility Xd/x* vs - T for a single hole in an spin background with 
J/t = 0.4 (full line), and J = (dashed line). 

and hence computational requirements are enhanced. On the other hand, we are evaluating only a 
thermodynamic quantity - free energy density T allowing for simplifications discussed in Sec. 4. 

We consider here only the case of a single hole iVj, = 1, doped into a magnetic insulator. 
Nevertheless we expect that results remain relevant for low Ch > 0. Namely, by assuming the 
independence of charge carriers (spin polarons) Xd should scale linearly with the doping, i.e. Xd °^ c h 
at least for Ch <C 1. We present in Fig. |6.10| results for Xd(T)/x* , as obtained by the analysis of a 
single hole in an AFM with J/t = 0.4 on a system with N = 20 sites. Here, x* = e^ta^jf? is a 
characteristic diamagnetic susceptibility, which can be e.g. compared with the Pauli susceptibility 
of tight-binding free fermions (with constant average density of states) XPi X* ~ ^Xp( m o/ m t) 2 
(Veberic et al. 1998). We calculate Xd(T) using the equation (|6.2l| ) and the difference between 
lowest fields, B = and B = Bq. For comparison we present also the J = case, where finite-size 
results can be checked with the HTE (Veberic et al. 1998). 

The J = case seems easier to study and to understand. The HTE and small systems show 
a gradual transition from the high-T regime of an incoherent hopping to the low-T Nagaoka FM 
state, accompanied with a monotonous increase of At T — > Xd is expected to diverge since 
within the g.s. the dependence on B is nonanalytic, i.e. 5E oc \B\ following from a cyclotron 
motion of a QP (FM polaron) in a finite B > 0. Still the asymptotic behaviour at low T is not 
simple to establish, since low lying states (above the Nagaoka state) are not easy to describe. 

The behaviour of the AFM polaron with J > is more involved. At T 3> J the exchange scale 
J is not important and results qualitatively follow that of J = 0. Quite remarkable is however the 
vanishing of the diamagnetic Xd (° r even the change of its sign) which appears on lowering T < t. 
This seems to follow from the fact that J > diminishes and even destroys emerging coherence of 
a QP. Only at lower T < J the coherence is reestablished and the standard dynamical picture of a 
coherent AFM polaron is dominating the behaviour in the lowest fields. Reliable results are however 
quite difficult to obtain since small systems results (note that we have to follow the variation with 
B) are quite sensitive to the system shape and boundary conditions due to a very anisotropic and 
degenerate QP dispersion. 
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At low Ch the relevant value of the diamagnetic susceptibility ChXd is small compared to the spin 
contribution (Veberic et al. 1998), hence it is not evident that it is directly measurable (Walstedt 
et al. 1992, Miljak et al. 1993). Nevertheless, the variation Xd(T) is in any case very interesting in 
connection with the challenging problem of the Hall effect, related via the equation ( |6.19| ), and the 
anomalous Rh(T). Both \d as we U as the Hall effect emerge from a coupling to orbital currents, 
which is the essence of the relation ( |6.19 ). We note that at T > J the HTE for the Hall constant 
Rh(T) is analogous to that of Xd(T) discussed here. Crossing the scale T ~ J remains the challenge, 
whereby it seems that at such T hole-like Rh > is even reduced from its high-T value (Shastry 
et al. 1993). Nevertheless Hall effect measurements (Ong 1990) indicate that Rh > recovers for 
T <T* , varying strongly with T, and finally seems to approach the well known quasiclassical result 
for T -» (Prelovsek 1997b). 



7 Spectral Properties 

One of the most desired quantities, giving valuable information on electronic properties in interact- 
ing systems, is the single-electron Green's function G(k,u) and the corresponding spectral function 
A(k,uj). In theoretical treatments G(k,u) is usually the basis to understand two-particle proper- 
ties, such as dynamical electric and magnetic response functions etc. In strongly correlated systems 
it appears however that calculations as well as the measurements of spectral functions are more 
delicate than that of most response functions, hence we also treat them towards the end. 

Spectral functions A(k, u) are directly accessible via ARPES experiments. There has been in 
recent years an impressive advance in the resolution of ARPES as well in the novel information 
on cuprates, mostly obtained on convenient layered BISCCO materials (see review by Shen and 
Dessau 1995, and references therein). In the regime of optimum-doped and overdoped cuprates 
ARPES experiments reveal for a wide class of materials a well defined large FS consistent with the 
Luttinger theorem and a QP dispersion close to a tight-binding band (Shen and Dessau 1995, Ding 
et al. 1996, Marshall et al. 1996). This seems to imply the validity of the concept of a metal with an 
electronic- like FS, although such a simple picture is in an apparent contradiction with anomalous 
magnetic and transport properties. The LFL interpretation is also spoiled by the overdamped 
character of QP peaks. Although a large background makes fits of particular lineshapes nonunique, 
the QP inverse lifetime is found to be of the order of the QP energy, supporting the possibility of 
the MFL scenario. 

On the other hand, in underdoped materials (Marshall et al. 1996) well defined QP crossing the 
FS gradually lose their character and spectral functions approach a qualitatively different regime of 
an undoped AFM (Wells et al. 1995), where the QP dispersion is that of a spin polaron in an AFM 
background. An excitement has been induced also by the leading-edge shift (Marshall et al. 1996), 
observed by ARPES in underdoped cuprates in the normal state. It indicates that a pseudogap 
consistent with the d— wave SC symmetry persists well above T c . On the other hand recent angle- 
integrated PES measurement (Ino et al. 1997b) indicate another higher energy pseudogap scale 
being in a closer correspondence with T*(qJ as deduced from xo{T), p(T) and Rh(T) (Batlogg et 
al. 1997). 

It has been unclear whether above features of spectral functions could be reproduced within 
generic models, such as the Hubbard and the t-J model, in particular in the most challenging regime 
of the intermediate doping. A(k,uj) in 2D models have been so far studied mainly via numerical 
techniques (Dagotto 1994), in particular using the ED (Stephan and Horsch 1991, Eder et al. 1994, 
Moreo et al. 1995) and the QMC (Bulut et al. 1994, Preuss et al. 1995, 1997). These studies have 
established a reasonable consistency of the model QP dispersions with the experimental ones, as 
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well as the possibility of a large FS, but have not been able to investigate closer the character of 
QP, as e.g. manifested in the corresponding self energies T,(k,cu) being in the core of anomalous 
low-energy properties. 

7.1 Green's function 

When applying the usual definition of the retarded Green's function (Mahan 1990) to the model 
of strongly correlated electrons, one should keep in mind that the t-J model ( |2.6[) acts within the 
restricted fermionic basis. To avoid the ambiguity it is convenient to build restrictions directly in 
the definition of retarded Green's functions using the projected operators ct , C£ s , being Fourier 

transforms of cj , Cj S , respectively, 



G R (k,uj) = -i I dt e - +t <{c^),4 s (0)} + >. (7.1) 



o 



Note that within the subspace of allowed states the definition is identical to the usual one. A(k, uj) 
can be represented as a sum of the electron spectral function A + (k, uj) and the hole spectral function 
A~(k,Lo), 

A(k,io) = —1mG R (k,u) = A + (k,oj)+A-(k,-uj), (7.2) 



7T 

which are expressed in terms of eigenstates as 



A + (k,uj) = n- 1 Y,\^ m \c\^n)\ 2 e-^ E --» N " ) 5(u + v-E m + E n ), 

n,m 

A-(k,-u) = n~ 1 Y,\(^n l \c l :^ n )\ 2 e~^-^8(-u J - fi - E m + E n ). (7.3) 



n,m 



A~ represents the dynamical response when one electron is taken out (hole added) from a system. 
It corresponds (within the sudden approximation) to PES experiments and is proportional to 
corresponding cross sections in ARPES. On the other hand, A + describes the dynamics on adding 
an electron and is related to IPES experiments. In the equilibrium electron and hole spectra are 
related via the Fermi function f{uS) = l/(exp(/?u;) + 1), e.g. A~(k, —uj) = f{uj)A{k,u). 

The omission of doubly occupied states has important implications, in particular we note the 
change in anticommutation relations, 

i i 

In the paramagnetic phase where (nj s ) = (1 — c/J/2, the equation ( |7.4j ) leads to a modified sum 
rule for the spectral function (usually normalized to unity), i.e. 

f°° -. 1 

a= duo A(k,u) = -(1 + c h ) < 1. (7.5) 



An important consequence of the relation ( |7.5| ) is an upper bound for the momentum distribution 
function, 

poo 

H = ( c h c ks) = / A~(k,uj)duj < a. (7.6) 
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The anomalous sum rule (|7.5| ) leads to some ambiguity in the definition of the self energy T,(k,uj). 
We keep the usual one, 

G R (k,u) = l/(w- (7.7) 

in order to retain the standard definition of QP parameters. An alternative would be to replace 
the nominator with a (Prelovsek 1997a). The equation ( [7.7] ) implies that the self energy does not 
vanish for \u\ — ► oo, but shows rather an asymptotic form T,(k,LO — > oo) — ► (1 — l/a)u. Note also 
that the t-J model doesn't allow for a free propagation even at J = 0, therefore we do not include 
any free term eg in the definition (|7.7|) . 

It is straightforward to calculate A{k,u) using the FTLM (Jaklic and Prelovsek 1997), again 
following the equation ( |3.42 ). Since operators ct ,cg s act between subspaces with different Nh, 



i.e. Nh — > iV/j =F 1, we in fact calculate separately A + and A~ according to definitions ( [7.3[) . In 
calculations at low T and fixed Ch = N^/N we also replace the grand-canonical average with the 
canonical one in the subspace of states with Nh holes. For a proper interpretation of results it 
is important to locate correctly the chemical potential /j,(T,Ch), discussed in Sec. 4.1. We have a 
nontrivial check of \x via the one-particle DOS 

M{e) = ^Y^A{k,e-n), (7.8) 

k 

which relates (j, to Ch, 

/oo 
M{u + ^){e^ + l)- l du = l-c h . (7.9) 
-oo 

We find a very good agreement between fi calculated from M(e) and the one determined from 
thermodynamics in Sec. 4.1. 

The advantage of the FTLM over the g.s. ED is that it gives smooth G R (k,Lo) even for quite 
low T < J provided that T > Tf s . This allows for a meaningful calculation of £(A;,u>), which gives 
a new insight in the properties of QP. We are thus able to evaluate also QP parameters defined via 
H(k,u>) (Mahan 1990), in particular the QP energy Eg, the weight consistent with the equation 
(|2.5| ), and the damping IV, 

E % = E'(k,Eg), 
Zf = l-9T!{%,u3)/du}\ u= E 1i , (7.10) 

r s = z^"(k,E % )\. 

7.2 Single hole in the antiferromagnet 

We begin by considering the spectral properties of a single hole injected in an AFM. Since in the 
undoped model particle-number fluctuations are suppressed, the chemical potential is not a well 
defined quantity in this case and we use here unreduced energy e instead of ui = e — [a. 



On the left side of Fig. 7A the spectral function of a hole in an AFM is presented. Although 
T = J/2 is relatively high, implying only short-range AFM correlations, still a coherent QP peak 
is clearly observed at low energies for most available k. For k = (tt/2,tt/2) and k = (27r/3, 0) 
peaks coincide with the g.s. minima for systems with N = 16 and N = 18 sites, respectively. 
This indicates that even at T > the major contribution to the coherent spectral weight comes 
from transitions between the g.s. configurations with Nh = and Nh = 1. The coherent QP peak 
shows a dispersion on the energy scale comparable to J, consistent with the self-consistent Born 
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Figure 7.1: Spectral functions for the undoped AFM (left), and for a system with finite doping 
Ch = l/N (right), both obtained on lattices with N = 16, 18 sites. 

approximation (Schmitt-Rink et al. 1988, Kane et al. 1989), which yields a bandwidth W ~ 0.6 t 
for J/t = 0.3 (Martinez and Horsch 1991). Since no long-range AFM order is expected in small 
systems, the propagation of the hole appears to be determined by short-range AFM correlations. In 
addition to the QP band at e < Eq ~ 2t, we observe also high-energy features at e Eq, which are 
only weakly dependent on k and can be attributed to the incoherent hole propagation. High-energy 
peaks have been observed also in g.s. ED studies, but the structure was claimed to disappear at 
large (Poilblanc et al. 1993). Our study shows a nontrivial structure consistent at all considered 
N = 16 — 20, so we would rather attribute peaks to resonance (excited) states of the AFM spin 
polaron. 
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In Fig. 7.2 we present the T-dependence of the spectral function at k* = (jr/2,ir/2), which 
corresponds to the g.s. momentum of a hole in an AFM. We realize from Fig. |7.2| that there is a 
significant spectral change going from T ~ to T > J. Namely, in the high-T regime the QP peak 
becomes progressively broader and merges with the incoherent background, which at the same time 
loses the detailed structure (resonance peaks). This development is plausible since for T > J we 
are dealing essentially with the hole propagation in a random spin background well accounted for 
within the RPA (Brinkman and Rice 1970). 
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Figure 7.2: Spectral function A(k*,e) for the undoped AFM at several T. 

At T < J the QP energy Er and the corresponding weight Z^ , as determined by equations 
(|7.10| ), show only a weak T dependence and we get Z = Z^{T) ~ 0.16. We note that this value 
is obtained for the spectral function normalized according to the relation (|7.5| ) with a = 1/2. It is 
easy to see that in order to make a comparison with the usual definition of the QP weight of an 
AFM polaron (Kane et al. 1989, Martinez and Horsch 1991) we should rather take Z = Z jo. and 
our value Z ~ 0.32 is consistent with the one Z ~ 0.284, obtained within the self-consistent Born 
approximation for J/t = 0.3. 

7.3 Spectral functions at finite doping 

In Figs. |7.1| , |7.3| we present the gradual development of spectral functions towards the optimum 
doping. For Ch ~ 0.06 and Ch ~ 0.12 results obtained on systems with N = 16, 18 sites are 
combined for = 1 and Nh = 2 holes, respectively. As distinct to the undoped case, the spectra 
are broadened with Lorentzians of variable width according to 5 = 5q + (Soo — <^o) tanh 2 (o;/A), with 
^oo = 0.2i, 5q = 0.04i, and A = l.Ot. In this way sharper (well resolved) lo ~ features remain 
unaffected, while some of finite-size structures at higher |u;| > t are smoothened out. In any case, 
5 is always taken smaller than the energy widths of main spectral features. 

At finite doping we observe in A(k,tu) at all available k a coexistence of sharper features, 
associated with coherent QP peaks, with a broad incoherent background, as already established 
in ED studies (Stephan and Horsch 1991, Moreo et al. 1995). Coherent peaks disperse through 
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Figure 7.3: A(k, uj) for = 2 holes on JV= 16, 18 sites, and for = 3 holes on N = 16 sites. 

uj = as k crosses the FS. Within the given resolution in the k space the FS appears to be large 
for all Ch > 0.12, consistent with the Luttinger theorem (Luttinger 1960). An analogous dispersion 
is observed for the underdoped case ~ 0.06, although the latter one requires a more careful 
interpretation as given lateron. The total QP dispersion W is broadening as Ch is increasing, 
qualitatively consistent with the slave boson result where W oc c^t + \ J (Baskaran et al. 1987, 
Wang et al. 1991). 

7.3.1 Intermediate doping 

We first discuss in more detail the regime of intermediate doping (Jaklic and Prelovsek 1997), in 



our lattices Ch ~ 0.12 and Ch = 3/16 ~ 0.19. In Fig. 7.4 we show T,(k,co) evaluated for Ch = 3/16 
at T = 0.1 t ~ Tf s . We first notice an asymmetry between the PES (uj < 0) and the IPES (uj > 0) 
spectra at all k. T,"(k,u) are small for uj > 0, as compared to uj < 0. For k outside the FS 
this implies a modest QP damping, consistent with sharp IPES peaks seen in A(k, uj) in Fig. |7.3|, 



containing the major part of the spectral weight. T,'(k,uj) shows an analogous asymmetry, in the 
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Figure 7.4: Self energy £(fc,u;) for Ch = 3/16 at various k. 

region uj > resembling moderately renormalized QP. Due to the definition ( |7.7|) the slope of X' is 
not zero even at |cj| S> t. 

The behaviour on the PES {uj < 0) side is very different. For all k, — are very large at 
<jj < —t, leading to overdamped QP structures. We should however distinguish two cases. For k 
well outside the FS, large |S"(o; < 0)| > t does not invalidate a well defined QP at uj > 0, but 



rather induces a weak reflection (shadow) of the IPES peak at uj < 0, as well seen in Fig. 7J3 for 
k = (vr,7r). On the other hand, for k inside or near the FS the variation with u is more regular, 
and can be directly related to the QP damping. Particularly remarkable feature in Fig. is a 
linear u dependence of T,"(oj < 0) for k = (vr/2, 0) and k = (it/2,tt/2). Meanwhile k = (0, 0), being 
further away from the FS, seems to follow a different more LFL-type behaviour. Similar general 
conclusions follow also from results obtained for the lower doping Ch = 2/18. 



To address the latter point in more detail, we show in Fig. 7.5 the T variation of T,"(k,uj) for 
both dopings at selected k inside the FS. For = 3/16 the linearity of S"(cj) is seen in a broad 
range —2t < u) < at the lowest T shown. Moreover, for this optimum doping the T dependence 
is close to a linear one, taking into account a small residual (finite-size) damping r/Q > at u) = 0. 
Data can be well described by 

-Y,"(k<k F ,Lo) = i lQ + 1 {\Lo\+tT), 7-1.4, £~3.5. (7.11) 



Such a dependence is consistent with one of the proposed forms (2.4) within the MFL scenario, 
as well as with the conductivity relaxation 1/t(oj) in the equation ( |5,17| ). In contrast, the T 
dependence for = 2/18 seems somewhat different, and we find — S" oc u only in the interval 
— t ^ w i$ — T. This would indicate the consistency with the original MFL form ( |2.3| ) (Varma et 
al. 1989), however we should be aware that in the underdoped regime finite-size effects are larger 
at given T. 
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Figure 7.5: E"(k,uj) for Ch = 2/18 and Ch = 3/16 at different T and selected k inside the FS. 



Here we should comment on the manifestation of the FS in small fermionic systems. At T, uj ~ 
we are dealing in the evaluation of the expression ( |7.3|) with the transition between g.s. of systems 
with Nh and N' h = ± 1 holes, respectively. Since these g.s. states have definite momenta &q, 
they induce sharp QP peaks for particular k = k' Q — kg defining in this way for a small system the 
FS, apparently satisfying the Luttinger theorem. At the same time the corresponding QP damping 
vanishes, i.e. Y,"(k,u ~ 0) ~ 0. 

From Y<(k,Lj) we can calculate QP parameters, as defined in equations ( 7.10 ). Results for 
Ch ~ 0.12 and Ch ~ 0.19 are presented by Jaklic and Prelovsek (1997). We note that parameters 
are of a limited meaning for k well inside FS due to a large damping T^. Otherwise, Er shows the 
enhancement of the dispersion with Ch, while < for > hp. To establish the relation with 
the FL theory one has to evaluate QP parameters at the FS k ~ kp. Of particular importance 
is the renormalization factor Z = (T = 0). We find that at lowest T ~ Tf s Zj: (T) is still 
decreasing as T is lowered, e.g. a variation (cca. 20%) is found within the interval 0.1 < T/t < 0.3. 
This is not inconsistent with the MFL form Z^ 1 ~ ln(u; c /T), although such a conclusion should be 
taken with care due to narrow T interval. 

Regarding the size of Z at low but finite T > we note, that the value of the momentum 
distribution function is very close to the maximum (7^) for all k < fcjr, i.e. ~ (1+c^) /2 (Stephan 



and Horsch 1991). Taking the FS volume according to the Luttinger theorem Vfs/Vbz = (1 — c ft)/2 
and assuming that falls monotonously with k, this implies an inequality for the discontinuity Z, 



5nr < 



2c h 



(7.12) 



kF 1 + Ch' 

One should keep in mind that a discontinuity is meaningful at T ~ 0, while at T > it indicates 
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only a gradual step. Still we find a consistent result Z = 0.28 for Ch = 3/16, while for Ch = 2/18 
the value Z = 0.35 is somewhat larger, possibly due to higher T used in calculations. 

An analogous argument can be used to explain the electron-hole asymmetry of A(k,u). Holes 
added to the system at k < kp, u < move in an extremely correlated system, strongly coupled to 
the spin dynamics (Prelovsek 1997a), the latter also exhibiting the anomalous MFL-type behaviour, 
as given by the expression ( |6.13| ). On the other hand, states for k > kp are not fully populated, 
allowing for a moderately damped motion of added electrons with to > 0. 

Let us comment here on the relevance of our results to ARPES spectra near the optimum 
doping. The main observation is that for k < kp we find the linewidth typically r ~ t ~ 0.4 eV 
well compatible with experiments at k away from the FS and at intermediate doping (Shen and 
Dessau 1995, Ding et al. 1996, Marshall et al. 1996). Also the MFL form has been claimed to 
describe better the experiments (Olson et al. 1990), although this point is still controversial (see 
Shen and Dessau 1995). 



7.3.2 Overdoped and underdoped regime 

Let us turn from the intermediate doping first to the overdoped regime. Spectra at higher doping 



Ch = 4/16 and Ch = 5/16 are shown in Fig. 7.6. While A(k,co) at = 4/16 are still qualitatively 
similar to the spectra at intermediate doping, at Ch = 5/16 they are already showing a substantially 
different behaviour. In the latter case the incoherent background in the PES spectrum is reduced 
while QP peaks for k < kp are sharpened and all of them are essentially underdamped with 
< E?. Nevertheless, we still have ~ E^ so it is not evident whether this is already a (more) 
normal LFL. 

The proper analysis of the underdoped regime is even more delicate. At Ch ~ 0.06 — 0.12, as 



presented in Figs. 7A, |7.3| , a 'shadow' feature is clearly seen in A(k,uj) for k > kp. Namely, along 
with the principal peak at to > a weak bump in the uo < part of the spectrum appears. In 
T,'(k,uj) for k > kp this effect emerges as a strong oscillation, most pronounced for k = (vr,7r). It 



leads even to a double solution of the equation ( [7.10 ) for E^ , analogous to the phenomena within 



the spin-bag scenario (Kampf and Schrieffer 1990). This effect is becoming less visible at larger 
doping, e.g. at Ch = 3/16, so its disappearance is related to the reduction of the AFM correlation 
length £ ~ 1. 

Let us try to make the contact with recent ARPES measurements on underdoped BISCCO 
(Marshall et al. 1996), which show the opening of the pseudogap even for T > T c , in particular 
near the k** = (n, 0) point, where the FS seems to disappear. On the other hand the QP peak 
disperses through the FS along the T-M direction where the FS seems to be well defined. If we 



look closer at Fig. 7.1 for the lowest finite doping ~ 0.06, we can recognize very similar features. 



We realize that close to k**, the PES part shows spectra very analogous to the undoped AFM 



in Fig. |7, 1 , where due to the doubling of the Brillouin zone the dispersion is expected to reach 
maximum at k* = (n/2,ir/2) (note however that on a 4 x 4 lattice with n.n. hopping k* and k** 
are equivalent) and folds back for k outside the reduced AFM zone. E.g., in Fig. [O] for Ch ~ 0.06 
we can clearly recognize PES peaks also for k = (ir,n/3) and k = (tt,tt/2), being lower in u than 
the QP peak at k*. Such a dispersion is also seen in ARPES. At the same time there is for both 
mentioned k > k* also a visible peak in the IPES part, which can be interpreted as a QP dispersing 
through the FS. The PES and the IPES peak in this case are separated by a gap, which does not 
seem to be a finite-size effect, and is qualitatively consistent with the ARPES feature (Marshall et 
al. 1996). It should be also noted that along the T-M direction the QP dispersion is closer to a 
normal metallic one with much less pronounced 'shadows' for k > kp. 
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Figure 7.6: A(&,u;) for Ch = 4/16 and = 5/16. 



7.3.3 Influence of next-nearest-neighbour hopping 



It is evident from spectral functions at the intermediate doping, e.g. in Figs. 7.3, |7.6| that the 
shape of the FS obtained within the t-J model is closer to a circular one than to the one found in 
cuprates via ARPES experiments (Shen and Dessau 1995), where e.g. k** is inside the FS. It has 
been well established that such a FS topology can be reproduced by adding the n.n.n. hopping term 
( |2.7| ) with t' < (Tohyama and Maekawa 1994). It is however still an open question whether the 
modified hopping and different FS topology are essential for the anomalous properties of cuprates. 



In order to see the qualitative changes introduced by the n.n.n. hopping, we present in Fig. 7.7 
A(k,uj) at fixed doping c/j = 3/16 for two different t'/t = —0.2 and t'/t = —0.35. It is evident 
that t' < lifts the degeneracy between of the k* and k** even in the 4x4 lattice. Also the QP 
dispersion changes and the QP peak with k* is now at ui > 0, i.e. outside the FS, while the QP 
with k** moves inside the FS. The effect is enhanced for larger \t'/t\. It is however important to 
realize that other QP properties are not essentially changed, in particular we find that the self 
energy H"(k,uS) remains qualitatively similar to the case with t' = 0. 
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Figure 7.7: A(k,u>) for c/j = 3/16 with nonzero n.n.n. hopping: t'/t = —0.2 and t'/t = —0.35. 



7.4 Density of states 



Finally we present in Fig. |7.8| the doping dependence of the single-particle DOS M{e) (Jaklic and 
Prelovsek 1997), as defined in the equation (|7.8|) . In a weakly doped system with Ch ~ 0.06, a QP 
coherent peak (of width ~ 2J) is seen at u> = e — fj, < 0. Besides that, a broad background due 
to well understood incoherent hole motion is dominating u <C 0. At such low doping the electron 
lo > part of the DOS is weak, with the total intensity 2ch as compared to 1 — Ch of the hole part. 
With increasing Ch the hole incoherent background does not reduce appreciably in intensity, while 
the coherent dispersion near the Fermi energy widens and cannot be well distinguished in M(e) 
from the background. At the same time, the electron part of the DOS is increasing, both in the 
weight and in the width. Note that oscillations which appear for uj > at higher doping Ch > 3/16 
are essentially finite-size effects. Namely we are dealing with a restricted number of finite-size k 
while QP peaks are becoming quite narrow. Such effects are even more pronounced in overdoped 
systems with c/j > 0.25, where the incoherent part is finally losing its intensity and the coherent 
QP are dominating the whole e regime. It should be also mentioned that the introduction of the 
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Figure 7.8: Single-particle DOS Af(e) for different dopings at fixed T/t = 0.15. For Ch ~ 0.06 and 
Ch ~ 0.12 we present joined DOS for N = 16, 18 systems with Nh = 1 and Nh = 2, respectively. 
The thin vertical lines denote the chemical potential fJ.(T). 



n.n.n. hopping t' ^ 0, in spite of a significant change of the FS shape, does not seem to induce an 
essential change within the DOS. 

Let us comment here the relation with the entropy s, Fig. |4.3| a. If we assume the low-T form 
as follows from the LFL theory (Abrikosov et al. 1965), we get 

s = T. (7.13) 

3Z 



We see from Fig. 7.8 that Af(p) is only weakly doping dependent at intermediate Ch and actually also 
quite close to the free-fermion value (assuming a tight-binding model with the hopping parameter 
t). Taking Z ~ 0.28 for Ch = 3/16, we get s ~ 0.29 ks at T = 0.1 t, quite consistent with 
thermodynamic results in Sec. 4.2. Here the surprising fact is that such s represents a large increase 
over the undoped AFM and appears due to a very low concentration of mobile holes introduced 
into an AFM. 

Let us now look closer on JV(e) in the underdoped regime. We realize quite clearly that at 
Ch = 0.12 a pseudogap starts to appear in the DOS &t u = e — /x~0 and becomes even more 



pronounced at Ch ~ 0.06. It is evidently related to a similar phenomenon in A(k,oj) in Fig. 7.1 
In order to establish the origin of such a gap, we follow in Fig. [7l] the variation of Af(s) with T. 
We find that the gap structure disappears at T > J, so it must be related to the onset of the 
short-range AFM order. Also the gap energy Ae < t seems to be determined rather by J. It is 
plausible that the onset temperature is related to pseudogap scale T*, discussed before in relation 
with the maximum in XoC^)j an d with the crossover in p(T). 

Recently the DOS has been measured via the integrated PES and IPES for LSCO compound in 
a wide range of doping < Ch < 0.3 (Ino et al. 1997b). More reliable are PES spectra proportional 
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Figure 7.9: jV(e) at various T for dopings: (a) = 1/20 and (b) Ch = 2/18. 



to J\f + (uj) = N{uj + /j,)f(uj), giving information on the hole part u> < 0. Some features are well 
consistent with our results. In the overdoped regime > 0.2 M + (uj < 0) appears quite flat. At 
the same time, there is a substantial and broad IPES contribution indicating an enhanced DOS at 
uj > 0. On the other hand, for Ch < 0.17 a pseudogap starts to emerge gradually at ui < 0. E.g., the 
inflection point in N + (uj) moves from uj ~ toward ui ~ —0.2 eV at lowest dopings. This is quite 
close to values found in Figs. Thus the gap scale indeed seems to be related to the pseudogap 
scale T* discussed in connection with \o- Note that the latter is substantially larger than the 
leading-edge shift ~ 30 meV found in ARPES (Marshall et al. 1996). Still it is well possible that 
both phenomena are closely related. 



8 Other properties 

8.1 Electron density correlations 

It is interesting to study also electron-density dynamics, as contained within the dynamical charge 
susceptibility Xc(q,^>) and the corresponding density correlation function N(q,u>), defined by 

MM = 4(l-e~^)N(M, 

POD 

N(q,uj) = Re / dt e^*(n f (t)n_ 9 -(0)), (8.1) 
Jo 

i 

Recently several numerical and analytical studies within the t-J model have been devoted to static 
N(q) (Putikka et al. 1994) as well as to dynamical N(q,u>) (Tohyama et al. 1995, Eder et al. 
1995, Khaliullin and Horsch 1996) at T = 0, mainly in connection with the interesting conjecture 
of the charge-spin separation (Anderson 1987, Baskaran et al. 1987) in layered cuprates. A parallel 
study of charge correlations N(q, ui) and the spin structure factor S(q, uj) however reveals quite an 
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essential difference between an ID Luttinger liquid and a 2D doped AFM (Tohyama et al. 1995), 
as well as between charge and spin fluctuations. Typically, N(q,u) show a broad peak at high 
a;~6tat<f~Q due to the incoherent motion of holes, while at q — > density fluctuations narrow 
in a collective acoustic-like mode (Khaliullin and Horsch 1996). In contrast to S(q,oj), N(q,oj) 
spectra seem to scale linearly with c/j at low doping (Eder et al. 1995). It should be also noted that 
quite generally N(q — * 0,u) fluctuations can be related to (anomalous) current correlations C(tu) 
through the particle-number conservation, so an anomalous behaviour is quite possible for N(q,oj) 
as well. 

The investigation of N(q, u) at T > has not been performed so far. We present here some 
results, using the FTLM in analogy with the evaluation of S(q,u) in Sec. 6. Our aim is again 
twofold. On one hand, FTLM gives smoother spectra at T > Tf s , which could be compared to g.s. 
ED results. More interesting is the low-frequency regime u < J, which cannot be studied reliably 
in g.s. calculations, and moreover the T dependence of spectra in this regime. 

In Fig. 8.1 we present normalized spectra N(q,ui)/ch for nonequivalent q on a lattice of N = 18 
sites and fixed T = 0.2 t. In this case we use J/t = 0.4. Shown are results for Ch = 1/18 and 
c/i = 3/18, belonging to regimes of low doping and optimum doping, respectively. We notice that 
normalized spectra are very similar for both Ch, at least for larger q, taking into account also that a 
more detailed structure for = 1/18 is partly due to finite size effects. Our results are in general 
close to previous numerical (Tohyama et al. 1995, Eder et al. 1995) and analytical (Khaliullin and 
Horsch 1996) results. It is characteristic that at q ~ Q the intensity is nearly exhausted by a large 
peak at u ~ 6t. As expected from conservation laws spectra are narrowing for q — > 0. 

It is however remarkable that in spite of quite low T in Fig. [O] N(q,uj — > 0) does not seem to 
approach zero for any q, the effect being more visible for smaller q. We follow the evolution of density 
fluctuations with T at intermediate doping Ch = 3/16 for chosen q = (7r/3,7r/3) in Fig. |8]^. It is 
evident that N(q,u) is nearly T independent even for u <T. The behaviour is quite similar to the 
one observed for local spin correlations S(u) in Fig. 6.5 and is consistent with the MFL assumption 



(|2.2|) for the charge susceptibility ^(u), requiring via the equation ( |8.1|) N(q,u> — > 0) = const 
independent of T. The MFL behaviour of N(q,iv) at intermediate doping is not entirely surprising 
due to the connection to the anomalous MFL-type C(uj) and optical conductivity ct(oj), as well due 
to the relation of the electron dynamics to anomalous spin fluctuations S(q,oj) (Prelovsek 1997a). 

N(q,u) have been mainly the subject of theoretical considerations so far. It should be however 
noted that a relation could be established with experimentally relevant dielectric function e(q, u), 
if the model would incorporate also the long-range Coulomb repulsion between electrons. Such 
an additional interaction could be possibly treated within a framework of an random-phase-like 
analysis. It would be desirable since e(q,u>) can be measured also in cuprates via the electron- 
energy loss spectroscopy (Niicker et al 1989) which yields directly Im[l/e(<f, to)]. 



8.2 Electronic Raman response 

One of very useful probes for the investigation of excitations in an AFM has been the Raman scat- 
tering. Using the latter method, it has been also established that the reference insulating cuprates 
correspond well to a planar Heisenberg AFM, where pronounced Raman resonance processes at low 
T are attributed to two-magnon excitations (Lyons et al. 1988, Singh et al. 1989). 

A general framework for the theoretical explanation of Raman processes in correlated systems 
has been so far given within the Hubbard model, where the effective Raman operator for resonant 
and off-resonant conditions has been derived (Shastry and Shraiman 1990). Still there are several 
unresolved questions concerning Raman processes in cuprates, and more generally in doped AFM. 
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Figure 8.1: Normalized density correlation spectra N(q,Lo)/ch at fixed T/t = 0.2 and J/t = 0.4 at 
two dopings: Ch = 1/18 and Ch = 3/18. 



One puzzling aspect concerns the observed pronounced T dependence of the two-magnon peak in 
undoped cuprates (Knoll et al. 1990). The latter has been interpreted as the phonon- induced 
broadening, but also invoking higher-order resonant processes (Chubukov and Frenkel 1995). An- 
other problem is the doping dependence of Raman spectra. Recent experiments, performed on 
YBCO materials in the resonant regime (Blumberg et al. 1994), show a dramatic increase of the 
broadening of the two-magnon peak with doping, so that spectra appear essentially flat in the 
normal phase when approaching the optimum doping. At the same time, the peak position does 
not move appreciably. 

Here we consider only the resonant Raman processes, since most experiments are performed in 
the resonant regime. In contrast to response functions considered in previous sections, the operator 
relevant for the resonant Raman scattering cannot be determined uniquely within the t-J model, 
since it necessarily involves higher (resonant) levels. Still we can adopt a view that more complete 
(e.g. three-band) model for cuprates can be mapped onto an effective Hubbard model (Hybertsen 
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Figure 8.2: N(q,u) for q = (ir/3, 7r/3) and different T at fixed c/,, = 3/18. 

et al. 1990) for resonant transitions of interest. Within the Hubbard model near half-filling the 
Raman-scattering operator has been derived in the limit t/U <C 1 (Shastry and Shraiman 1990) 
yielding the well known form for the Heisenberg AFM (Parkinson 1969), 

R = A^2(e inc ■ ?ij)(e sc ■ fij)(Si ■ Sj - jr^n,), (8.2) 

where e*j nc , e sc are the incident and the scattered electric-field unit vectors, respectively, and the 
amplitude factor A = At 2 / (U — uii nc ) oc J incorporates the resonance at the incident-light frequency 
^inc ~ U. The Raman spectral function is then given by 

I(u) = Re / dt e lult (R(t)R(0)}, (8.3) 

which can be calculated at T > using the FTLM analogous to other correlation functions. 

Within the t-J model with finite doping the Raman intensity I(oj), corresponding to the operator 
( |8.2| ), has been evaluated so far at T = using the ED (Dagotto and Poilblanc 1990), while for 
the undoped AFM also T > has been studied using the full diagonalization (Bacci and Gagliano 
1991). Results obtained via the FTLM have an advantage over previous ones also for low T, since 
the Raman intensity is not expected to be singular for T — > 0, so much smoother spectra are 
obtained by using small but finite T > Tf s . 

We restrict our analysis to the dominant B\ g scattering geometry with e*j nc = (e x + e y )/y/2 and 
e*sc = (ex — e y )/V2. The undoped AFM at T > Tj s ~ J/2 has been studied by Prelovsek and Jaklic 
(1996), and reveals at low T a two-magnon Raman peak at uj ~ 3.3 J, consistent with experiments. 
The width is however quite narrow and starts to broaden substantially only at higher T ~ J, where 
a gradual transition to a broad featureless spectrum occurs. Hence other mechanisms have to be 
invoked to account for the observed pronounced T-dependent width at lower T -C J (Knoll et al. 
1990). 
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In doped systems T plays a less essential role and typically we observe only a weak decrease 
of the Raman intensity in the interval T/J = 0.3 — 1.0 (Prelovsek and Jaklic 1996). On the other 
hand, the dependence on doping is essential, as evident from Fig. |8.3| , where we present I(u>) for 
various dopings Ch < 0.25 at lowest T = 0.15 i > TV S . Already the smallest nonzero doping 
c/j = 0.05 increases dramatically the width of the two-magnon peak, while spectral features become 
overdamped on approaching the optimum doping Ch ~ 0.15. It is however remarkable, that the 
peak position does not shift appreciably in the underdoped regime < 0.1. Only for > 0.15 
the spectra change to a broad central-peak form with a maximum at u) = 0. 
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Figure 8.3: Raman intensity I(uj) for different dopings c/j at fixed T = 0.15 t. 



How can we interpret the above results for the B\ g resonant Raman scattering ? Up to the 
overdoped situation ~ 0.3 the scattering seems to be dominated by the spin-exchange part flS.2|) , 
since the latter determines the low-energy fluctuations even at the optimum doping. Still, there 
are evident changes with doping, since the reduced AFM correlation length at larger induces a 
large broadening of the Raman two-magnon peak and a reduction of its intensity. In this respect 
the doped system behaves quite similarly to an undoped AFM but at elevated T ~ J. 

A systematic resonant-Raman scattering study has been recently performed on a sequence of 
YBCO materials (Blumberg et al. 1994). It seems that in the underdoped regime our model 
results well reproduce experimental ones, both regarding the shape of the Raman spectra and their 
intensity variation with doping. At optimum doping experiments still reveal a weakly pronounced 
peak at the same frequency u ~ 3 J, while our results on Fig. 8.3 already reveal a maximum at 
uj = 0, nevertheless both spectra are nearly flat. 



8.3 Thermoelectric power 

Among other transport properties the thermoelectric power (TEP) or the Seebeck coefficient S is 
one of the most frequently investigated. Within the cuprate family TEP has been measured for 
several materials, for an overview see e.g. Kaiser and Uher (1991) and more recent analyses (Tallon 
et al. 1995, Cooper and Loram 1996). Again an universal behaviour has emerged depending mainly 
on hole doping. For weakly doped materials S is large and positive, with some T dependence at 
higher T. With increasing doping TEP S decreases rapidly. It also falls off nearly linearly with 
T, but with a rather small slope. Quite consistently S is nearly vanishing for the optimum doping 
and changes sign to a negative one in overdoped materials. This variation has clearly some parallel 
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with the doping dependence of the Hall coefficient Rh, also showing a transition from a hole- like 
behaviour Rh > at low doping into an electron-like Rh < in overdoped samples. It seems 
plausible that both anomalous properties are related and emerge from strong correlations in these 
materials. 

Theoretically TEP in doped AFM or more generally in strongly correlated metals has not 
attracted much attention so far (Hildebrand et al. 1997). This is not surprising, since it seems to 
require the understanding of both electrical and thermal currents. 

Within the linear response theory (Mahan 1990) the TEP can be expressed in terms of particle 
current j and energy current je correlation functions, 



s 1 



e T 



M C(u -» 0) 



(8.4) 



where mixed correlation function Cj E j(uj) is defined in analogy with the expression (5.1) for C(uS). 

The evaluation of S(T) thus requires the calculation of Cj E j(u>), in addition to the current- 
current correlation function C(u) discussed extensively in Sec. 5. It is straightforward to derive 
the expression for jg, nevertheless the operator is more involved due to three-site terms. We do 
not intend to give here a more complete analysis of Cj E j(u>) which will be presented elsewhere. 
We mention however only a preliminary observation, that Cj E j(u) and C{oj) appear to be closely 
related, 

C 3e3 {uj) = (C(u), (8.5) 

i.e. ( is nearly u> independent for u < t, but as well T independent at low T < J. Since we are 
dealing at Ch > with a metal, although with a strange one, we expect a finite \S(T — > 0)| < oo. 
Hence we deduce from the equation fl8.4j ) the condition ( = fj,(T = 0). Our numerical results for £ 
are indeed consistent with this assumption. 
We thus arrive at the simplified expression, 

S~^[» h (T = 0)-v h (T)], (8.6) 

which involves only /^(T), studied in Sec. 4.1 and shown in Fig. |4.1|, In the low-T regime of interest 
we claim in a broad regime of doping < 0.3 a linear variation (|4.3| ) of /x^(T). This leads directly 
to S ~ — aSo where Soks/eo = 86/jlV/K. It is evident from Fig. |7.3| that a is changing sign from 
a negative for < c* h to a positive one for > c* h . In Fig. |8.4j we plot our result for S{ch). 
Note that due to simplifications involved in the expression ( |8.6| ) we do not intend to consider the 
T variation, hence our results apply to lowest reachable T < Tf s ~ 0.1 t. For comparison we 
present on the same plot also experimental results for LSCO and oxygen deficient YBCO, taken 
from Cooper and Loram (1996), whereby data refer to the normal state at T ~ 300-fT. Qualitative 
agreement between theory and experiment is quite reasonable. In LSCO values of S are decreasing 
with doping, while S is quite high at low doping Ch < 0.05. Towards optimum doping S essentially 
vanishes. Similar is the trend for YBCO data, although there are clearly quantitative differences. 
The main disagreement between the calculated and experimental S is in the overdoped regime, 
where our results indicate S < with probably too large values. It is however well possible that 
in this regime our analysis is not adequate, first due to ^h{T) approaching more normal LFL T 2 
dependence, as well as due to the breakdown of the relation ( |S.5| ). 
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Figure 8.4: Thermoelectric power S vs. Ch for T/t = 0.1. Experimental result for LSCO and 
oxygen deficient YBCO are taken from Cooper and Loram (1996). 

9 Discussion 

In the absence of accepted microscopic or phenomenological theories, which would at least quali- 
tatively describe normal-state properties of cuprates in the whole doping regime, it is also hard to 
summarize our numerical results in a compact manner. 

It is first important to understand which phenomena determine the optimum doping in cuprates 
and in microscopic models. In cuprates the optimum is defined by T c {c° h pt ) = max, but at least 
very close in doping, i.e., c\ ~ c^ 1 are also materials where the low-T electronic entropy s and the 
uniform susceptibility xo are maximum (Loram et al. 1996), where the pseudogap scale disappears, 
i.e. T* -> T c (Batlogg 1997), etc. 

Within the t-J model we can associate the optimum so far only with respect to normal-state 
properties, e.g. with -s(c^) = max. Such a maximum has to exist at a nontrivial c* h > due to 
the relation (|4.6| ) with /2h(T), where the latter must change the sign from a weakly doped AFM 
(semiconductor-like regime) at Ch > to a regime of nearly free 2D fermions at c e ~ 0. It is plausible 
that within the prototype model ( [2.6D c k is determined by the interplay of the spin exchange energy 
Ej = (Hj) and the (hole) kinetic energy E t = {H t ). Since at low doping we get \E t \ oc tch and 
\Ej\ oc J we can get a rough estimate c\ oc J/t <C 1. From such an argument it seems plausible 
that at Ch ~ c* h the AFM correlation length becomes very short £ ~ 1, while AFM correlations are 
essentially irrelevant in the overdoped regime > c* h . 

9.1 Universality at intermediate doping 

It is evident from our results that normal-state properties are universal at the optimum doping 
Ch ~ c* h . Since the maximum s(c/J is broad within the t-J model also this regime is quite extensive, 
i.e. 0.15 < Ch < 0.25 for chosen J/t = 0.3. The main feature of this regime obtained in our 
study is the MFL-type dynamics in several response functions: the local spin susceptibility Xl(<^)) 
the optical conductivity ct(lu), density correlations N(q,to), and spectral functions A(k,co). It is 
characteristic that the dynamics in the low-u; (lu < J) regime seems to be universal, i.e. it is 
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determined solely by T itself without any additional parameters, as concluded from equations 
( 5.131 ) for cr(u) and (|6,13j ) for x'l( u )- ^ should be however reminded that this universality could 



be restricted to a certain T range, e.g. we can show its validity only for Tf s < T < J. 

It seems that the origin of the universality is in the large degeneracy of the spin subsystem, 
where the frustration in induced by the hole doping. A plausible argument is that holes tend to 
induce FM correlations to minimize kinetic energy, while spins prefer an AFM ordering. In any 
case, only spins have enough degrees of freedom to explain large entropy s ~ 0.4 ks even at T < J. 

Let us attempt an explanation for the universality of S(lo) at ch ~ c* h . From the explicit 
representation in terms of eigenstates, 

S(u) = (1 + e"^) — ^— K*«l^ |*m)| 2 % - E m + E n ), (9.1) 

n,m 

we can first discuss conditions that the response is T independent for uj > T. Since in this case 
the prefactor is constant, a naive requirement appears to be that low lying states E n — E$ < T 
have similar matrix elements to excited states with E m — E n > T. Nevertheless the validity of this 
observation clearly depends on the character and on the density of low-lying many-body states, 
and e.g. does not apply at low T to a system with a unique g.s. As realized from the entropy 
results in Sec. 4.2 the degeneracy and the density of many-body states is clearly very large at the 
intermediate doping. To explain the T independence of S(u>) even for uj < T, we recall the sum rule 
flfTTiD and assume that there is no characteristic scale uj c < T which could introduce an additional 



low-w structure in S(u). A natural candidate for such a scale in an AFM is the (gap) frequency 
w c ~ c/£ where c is the spin wave velocity and £ the AFM correlation length. Here originates the 
essential difference between an undoped and the optimally doped AFM. In a pure AFM within the 
renormalized classical regime £ is exponentially large for T -C J (Chakravarty et al. 1989) and 
consequently uj c <C T. On the other hand, at the intermediate doping £ appears to be short and 
not strongly T dependent, i.e. £ oc l/^/cf t . The latter situation excludes u c < T < J, hence leading 
to an universal S(u). 

Turning to the discussion of a(oj), as given by the equation fl5.13|) , we realize that it is not 
meaningful to associate l/r(0) oc T with a current relaxation rate. Namely the current-current 
correlation function C(u>) shows a broad uniform spectrum indicating on a very fast current decay 
with 1/f ~ 2t. Thus it seems that a general condition for the validity of the specific MFL form 
( |5.13| ) is a fast current decay with 1/f 3> T. This could happen due to holes moving in a disordered 
spin subsystem, which leads to entirely incoherent collisions among holes and the appropriate mean 
free path l s is only few unit cells. Hence it seems that we are dealing with a novel phenomenon of 



quantum diffusion which follows an universal form ( 5.13 ). 

In an analogous way also the relaxation of single electron states near FS is affected, as manifested 
in the MFL form of the self energy H(k,uj), equation (|7.11|) . Anyhow it is plausible that the QP 
damping T,"(k ~ kp,uj) and the current relaxation 1/t(u), as given by the equation ( |5.17| ), are 
closely related. There are also phenomenological arguments (Littlewood and Varma 1991) which 
indicate that the MFL-type dynamics of a boson (spin) subsystem implies a MFL behaviour of the 
electron propagator. Recently such a relation has been derived by one of the authors (Prelovsek 
1997a), employing a decoupling approximation for the projected G R (k,u). 

On the other hand it should be stressed that within the same regime certain electronic properties 
appear quite normal, i.e. close to the LFL behaviour. First, spectral functions indicate on a rather 
well defined large FS following the Luttinger theorem. This seems to indicate that correlated 
electronic state evolves continuously from a noninteracting one, whereby the FS can be traced to 
momenta of lowest (g.s.) many-body states. It is much harder to explain quite LFL behaviour of 
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the entropy s oc T at low T, as seen in Fig. 4.3a, and moreover a free- fermion- like Wilson ratio 
W ~ Wo, as found in Sec. 6.2. 

It should be noted again that such an universality could be restricted to a certain T window. 
Namely it is well possible that at low T ~ T co h some coherence or ordering could appear. Note 
however that e.g. the onset of a spin ordering would invalidate our arguments on the MFL-type 
universality in the low-T phase. This is relevant for cuprates showing stripe structures (Tranquada 
et al. 1995) or longer range incommensurate ordering (Hayden et al. 1996). Even in this case the 
dynamics at higher-w is possibly not affected by the ordering, as seen in the analysis of the neutron 
scattering experiments in Sec. 6.6. 

At the intermediate doping our numerical results do not give an indication for a coherence 
temperature T co h, where the breakdown of the universality is expected. This is not surprising in 
view of the experimental facts for cuprates where at optimum doping ~ T c and our Tf s > T c . 
Nevertheless we note that certain restrictions follow already from the MFL-type dynamics. E.g., 



the particular MFL form (5.13) for cr(uj) leads at T — > to a divergent integral on the lower bound 
u) — » and hence violates the optical sum rule ( |5.S| ). Therefore such a form cannot apply to 
arbitrary low T giving an indication for the existence of a lower crossover T co h- 

9.2 Energy scales in underdoped antiferromagnet 

The underdoped regime Ch < c* h is less convenient for the FTLM approach. One reason is that 
Tf s is increasing towards the undoped AFM. At the same time the low-energy (pseudogap) scale is 
increasing and its effects are for certain quantities hardly resolvable from finite-size effects, hence 
they increase the uncertainty of results. 

Still there are clear indications for the pseudogap scale T* within the t-J model. Most evident is 
the maximum in Xo{T), as presented in Sec. 6.2. Quite a similar behaviour can be followed also for 
s(T)/T, as deduced e.g. from Fig. |4.3| a. This is related to nearly constant Wilson ratio Rw(T) ~ 1, 
as defined in the equation (|6.7D . Note that such a behaviour has been found also experimentally by 
Loram et al. (1996), extending to low-T regime not reachable in our calculations. It is quite evident 
from x(T) that the pseudogap scale T* is related to the onset of short range AFM correlations. 

It seems that the same phenomenon governs the pseudogap appearing in the DOS M(e) in the 
underdoped regime, as discussed in Sec. 7.4. While the pseudogap is clearly visible for lowest Ch, e.g. 
for Ch ~ 0.06, it is becomes more shallow towards the optimum doping ~ ct. This is consistent 
with PES measurements by Ino et al. (1997b). Experimentally one of the indications for T* is also 
the change of the slope in p(T). Here our results are less conclusive and appear to be finite-size 
limited for T < T*. A possible interpretation is that the QP scattering in the presence of longer 
range AFM correlations becomes a different one, in particular the mean free path l s increases due 
to the opening of the spin pseudogap. In our small systems this possibly leads to l s (T < T*) > L 
and consequently to a finite-size D C (T < T*) > which cannot be interpreted uniquely. Note that 
so far there is neither a satisfactory theory nor a numerical calculation of the mobility /j*o(T < J), 
even for a single hole in an AFM which could serve as a guide for the limit Ch — > 0. 



Experiments indicate the existence of a lower crossover scale T sg (Batlogg 1997), see Fig. 2.1 
identified as a spin gap in the NMR relaxation and appearing also in the ARPES, ct c (lj) etc. In our 
calculations it would be easiest to resolve the existence of a spin gap in Xo(T) and in S(u>), still we 
do not find an indication for that down to T ~ 0.1 t. 
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9.3 Conclusions and open questions 

One of the main conclusions of this work is that the prototype t-J model, in spite of its apparent 
simplicity, accounts surprisingly well for a variety of normal state properties of cuprates. This 
confirms the belief that most unusual properties are dominated by strong correlations and by the 
interplay of antiferromagnetism and itinerant character of electrons, both effects being inherent 
within the model. The agreement is a qualitative but as well a quantitative one, hence the param- 
eters of the model appear to be really microscopic ones in the whole range of doping and not just 
some effective parameters changing with doping. 

There are evidently many open questions regarding the interpretation and the understanding 
of our results, and related experimental facts on cuprates: 

a) The universal MFL-type dynamics in the intermediate regime, which seems to be well founded 
by our results on ct(uj), x'ii^) an d self energies, as well as by experimental facts, lacks a proper the- 
oretical explanation, although some relations have been proposed already (Littlewood and Varma 
1991, Prelovsek 1997a). 

b) There are more fundamental problems in the underdoped regime. One of them is the evolution 
of the FS in the limit Ch — * 0. A number of authors investigated the evidence for the small FS 
- hole pockets (Eder and Becker 1991, Eder et al. 1994) and its possible consequences (Trugman 
1990) at low doping, also stimulated by recent ARPES measurements (Marshall et al. 1996) of 
underdoped materials. In a strict sense at T — > the transformation from a large to a small FS 
seems to require a phase transition with a qualitative change in a number of properties. There is 
no evidence for that neither in experiments nor in our small-system analysis. On the other hand 
our spectral functions at lowest > reveal a folding of QP dispersion analogous to a single hole 
in an AFM. Still it is possible that such a feature coexists with a normal QP dispersion crossing of 
a large FS, however with a weight Z — > 0, i.e. gradually vanishing with doping as predicted by the 
equation ( 7.12| ). Such a scenario is not in contradiction with ARPES results. It has been in fact 



recently deduced from the integrated PES experiments (Ino et al. 1997b), still it would be hard to 
establish it beyond doubt. 

c) Related is the question of the thermodynamics at very low doping. There is a distinction whether 
at low T QP behave as a degenerate Fermi gas with finite QP weight Z > even for Ch — > 0, or as a 
nondegenerate system of spin polarons with Z — > 0. In the first case one could expect for Ch, T —* 
a finite electron compressibility k < oo, as discussed in Sec. 4.1, while the alternative would induce 
K — > oo. Our results are not conclusive in this respect, there is however also a numerical evidence 
for the latter scenario within the Hubbard model (Furukawa and Imada 1992, Assaad and Imada 
1996). 

There are several quantities which we did not discuss in this review. One of most challenging 
is the Hall effect, which is known to be anomalous. Unfortunately the inclusion of a real magnetic 
field B > increases computational efforts (see Sec. 6.7), while at the same time also the physics 
of the Hall conductivity is less local. The anomalous properties of the orbital diamagnetism Xd(T) 
and its relation to the Hall effect, equation ( |6.19 ), nevertheless confirm the anomalous behaviour 



of Rh(T). Also perpendicular a c {uS) is known to be a challenge for theoreticians. Still it appears 
that the interplane transport is essentially incoherent for T > T c , so it could be possibly related to 
the knowledge of planar spectral functions A(k,uj). 

The major open question within the t-J (or Hubbard) model is still the existence of additional 
low energy scales and of related low-T transitions. There are certain indications that such scales 
must exist in such models, it remains however a subject of future studies to find whether such 
transitions lead to an ordered spin and charge structure or to desired superconductivity. 
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